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Abstract 

In this paper, we consider TV-dimensional real Wishart matrices Y in the class 
Wr (£, M) in which all but one eigenvalues of £ is 1. Let the non-trivial eigenvalue 
of S be 1 + t, then as N, M — > oo, with M/N = ~/ 2 finite and non-zero, the eigen- 
value distribution of Y will converge into the Marchenko-Pastur distribution inside 
a bulk region. When r increases from zero, one starts to see a stray eigenvalue of 
Y outside of the support of the Marchenko-Pastur density. As the this stray eigen- 
value leaves the bulk region, a phase transition will occur in the largest eigenvalue 
distribution of the Wishart matrix. In this paper we will compute the asymptotics 
of the largest eigenvalue distribution when the phase transition occur. We will first 
establish the results that are valid for all N and M and will use them to carry out 
the asymptotic analysis. In particular, we have derived a contour integral formula 
for the Harish-Chandra Itzykson-Zuber integral Jq/^\ e tr ( XaY9 ^g T dg when X, Y are 
real symmetric and Y is a rank 1 matrix. This allows us to write down a Fredholm 
determinant formula for the largest eigenvalue distribution and analyze it using or- 
thogonal polynomial techniques. As a result, we obtain an integral formula for the 
largest eigenvalue distribution in the large N limit characterized by Painleve tran- 
scendents. The approach used in this paper is very different from a recent paper [23], 
in which the largest eigenvalue distribution was obtained using stochastic operator 
method. In particular, the Painleve formula for the largest eigenvalue distribution 
obtained in this paper is new. 



1 Introduction 

Let X be an N x M (throughout the paper, we will assume M > N and N is even) matrix 
such that each column of X is an independent, identical iV-variate random variable with 
normal distribution and zero mean. Let S be its covariance matrix, i.e. = E(XnXji). 
Then £ is an N x N positive definite symmetric matrix and we will denote its eigenvalues 
by 1 + Tj. The matrix Y defined by Y = j^XX T is a real Wishart matrix in the class 
Wr (£, M) and M is called the degree of freedom of the Wishart matrix. We can think of 
each column of X as a draw from a iV-variate random variable with the normal distribution 
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and zero mean, and Y the the sample covariance matrix for the samples represented by X. 
Real Wishart matrices are good models of sample covariance matrices in many situations 
and have applications in many areas such as finance [19], |JT], genetic studies and climate 
data. (See [IB] for example.) 

In many of these applications, one has to deal with data in which both N and M are 
large, while the ratio M/N is finite and non-zero. In this case, the sample covariance matrix 
Y becomes a poor approximation for the covariance matrix E. However, since it is often 
reasonable to approximate the data by iV-variate gaussian random variables, a comparison 
between the eigenvalue distribution of the sample covariance matrix and Wishart matrices 
with a known covariance matrix will give a good estimate for the spectrum of the true 
covariance matrix E. In particular, in applications to principle component analysis, one 
would like to study the asymptotic behavior of the largest eigenvalue of Y as N, M — > oo 
with M/N -> 7 2 > 1 fixed. 

For many statistical data with large N and M and M/N finite, it was noted in [IB] 
that the eigenvalue distribution of the sample covariance matrix is well-approximated by 
the Marchenko-Pastur law [5T] inside a bulk region (See [S],[B], [7], [8], [02] ■) 



where is the characteristic function for the interval and b± = (1 ± 7 -1 ) 2 . 

However, outside of the bulk region, there are often a finite number of large eigenvalues at 
isolated locations. This behavior prompted the introduction of spiked models in [18], which 
are Wishart matrices with a covariance matrix E such that only finitely many eigenvalues 
of E are different from one. These non-trivial eigenvalues in E will then be responsible for 
the spikes that appear in the eigenvalue distribution of the sample covariance matrix 
The number of these non-trivial eigenvalues in E is called the rank of the spiked model. 

Of particular interest is a phase transition that arises in the largest eigenvalue distri- 
butions when the first of these spikes starts leaving the bulk region. This phenomenon was 
first studied in [TD] for the complex Wishart spiked model and then in [7T] for the rank 1 
quarternionic Wishart spiked model. Let 1 + r be the non-trivial eigenvalue in E, in both 
cases, it was shown that the phase transition in the largest eigenvalue distribution occurs 
when r = / y~ 1 . Their results are expressed in terms of the Hastings-McLeod solution of 
the Painleve II equation, which is the unique solution to the Painleve II equation 



P(A) 



2ttA 



y/(X-b-)(b+- X)x[b.,b+], 



(1.1) 



0o(C) = C0o + 2$, 



(1.2) 



with the following asymptotic behavior 



0o ~ Ai((), C -> +oo 




(1.3) 



In particular, they have proved the following. 
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Theorem 1.1. (HUtf for complex and f7T$ for quarternions) Let W(E, M, 2) = II ' •(>]. M) 

and W^(E,M, 4) = Wq(E, M) be the N -dimensional complex and quaternionic Wishart 
matrices with M degrees of freedom respectively. Suppose all but one eigenvalues of E are 
1 and the other eigenvalue is 1 + r. Then as N , M — >■ oo with M/N — > 7 2 > 1 finite, we 
have 

1. For — 1 < r < 7 _1 , i/ie largest eigenvalue distribution is given by 

Jim P [ (A W(KC - (1 + 7- 1 ) 2 ) < C J = TW P (C) 

(1 + 7)3 / 



M-5-oo 



where TWp(Q is the Tracy-Widom distribution. 

TW 2 (() = exp (- jHfo - Q<l>l(y)dy 

TW 4 {() = -\/TW 2 {C,) (e~^ J c°° ^ + e \ ST 
2 V 



2. For r = 7 l , 



. 2 



lim p ( (A mai . - (1 + 7- 1 ) 2 ) < c ) = m) 

' ~ (l + 7) r 



M->oo 



where Fp(Q is given by 



F 2 (0 = iwf(C), F4(C)=t^(C), 7W X (0 = VwRCK 



2 / c °° <My 



5. For r > 7 _1 ; 

IT P f( A — - (^ + !) C 1 + T^" 1 )) 7: 2 2 <()= erf (C) 

A/^oo (1 _|_ r )yi _ T -2 r - 2 / 

where erf(() is the error function. 



The functions TW^g are called the Tracy-Widom distributions in the literature and 
they give the largest eigenvalue distributions for the real ((3 = 1), complex (/3 = 2) and 
quarternionic (/3 = 4) Wishart ensembles with E = I, as well as the largest eigenvalue 
distributions for a large class of random matrix models. 

Note that in [10], the phase transition was in fact computed for spiked models of any 
finite rank. These previous results naturally divides the range of the non-trivial eigenvalue 
1 + r into 3 regimes, which are illustrated in Figure [TJ 
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Figure 1: The 3 different regimes. Left: The subcritical regime where the perturbation in 
£ is not strong enough to form a spike in the spectrum of Y. Middle: The critical regime 
where the perturbation is just strong enough to form a spike at the edge of the spectrum. 
Right: The super-critical regime where a spike has left the spectrum 

1. The subcritical regime: When — 1 < r < 7 -1 , the perturbation in S is not strong 
enough to form a spike in the eigenvalue distribution of Y and the largest eigenvalue 
distribution in the Wishart matrix is not affected by this non-trivial eigenvalue. For 
real Wishart ensembles, this case was studied in [39] and it was shown that the largest 
eigenvalue distribution remains the same as the case when E = /. The result of [31] 
also applies to much more general sample covariance matrices that are not necessarily 
Gaussian. 

2. The critical regime: When 1 — 7~V = O ^M~3^, the perturbation is just strong 

enough to form a spike and a phase transition occurs in the largest eigenvalue dis- 
tribution. The result in JTU] for the complex phase transition in fact extends to the 
whole regime where the authors showed that the largest eigenvalue distribution is 
described by a function F 2 ((,w) which equals TWf(£) when w = 0. The critical 
regime for the real case is the subject of this paper. It was also studied recently in 
[23] using a completely different approach. 

3. The super-critical regime: When r > 1, the perturbation in S is strong enough to 
form a spike and the largest eigenvalue is located at the spike instead of the edge of 
the bulk region. For rank 1 real Wishart spiked model, this was studied in [55]. The 
results in [SB] has also been generalized to a large class of rank 1 spiked perturbed 
random matrix models in [72J. 

Despite having the most applications, the phase transition for real Wishart spiked model 
has not been solved until very recently [23] . The main goal of this paper is to obtain the 
asymptotic largest eigenvalue distribution for the rank 1 real Wishart spiked model in the 
critical regime. In a recent paper [23], the asymptotic largest eigenvalue distribution for 
the rank 1 real Wishart ensemble was obtained by using a completely different approach 
to ours. In [23], the authors first use the Housefolder algorithm to reduce a Wishart 
matrix into tridiagonal form. Such tridiagonal matrix is then treated as a discrete random 
Schrodinger operator and by taking an appropriate scaling limit, the authors obtained a 
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continuous random Schrodinger operator on the half-line. By doing so, the authors in [23] 
bypass the problem of determining the eigenvalue j.p.d.f. for the real Wishart ensemble 
and obtain the largest eigenvalue distribution in the asymptotic limit. In [23], the largest 
eigenvalue distribution is characterized in several different ways, one of which is the solution 
of a PDE. 

Theorem 1.2. (Theorem 1.7 of JEM) Let w = ( (1+^1)2 ) " i 1 ~ 7 r ) G (-°°> then the 
following boundary value problem 

OF 2d 2 F .. 2\ 9F 

F((,w) 1, (,w — > +oo together, 
F((, w) — > 0, w — > — oo, x < xo < oo 

has a unique solution Fp((,w) and 

hm P[ (A mag -(l + 7 - 1 ) 2 ) 7M3 <C) =F fi {C,w), 

(1 + 7)3 I 



M—>oo 



where (3 = 1,2 and 4 for the real, complex and quarternionic Wishart ensembles respec- 
tively. 

The distribution also has another characterization which requires more set up to explain. 
We will refer the readers to [23] for more details. The above result in [23] in fact also extends 
to the spiked perturbations of Gaussian ensembles in random matrix theory and to their 
general (3 analogues, which are defined in [23J. 

On the other hand, the approach in this paper uses orthogonal polynomial techniques 
that are closer to those in [ID] and [71] and we obtain a characterization of the largest 
eigenvalue distribution in the critical regime in terms of Painleve transcendents. It will 
be very interesting see how these two representations of the largest eigenvalue distribution 
can be converted into one another. We will now state our results and outline the method 
used in the paper. 



2 Statement of result 

We will now explain the approach used in this paper. 

Let Xj be the eigenvalues of the Wishart matrix, then the j.p.d.f. for the real Wishart 
ensemble is given by 

1 N - - r 

P(A) = -^|A(A)inAr^ / e-f^ Y ^g T dg, (2.1) 

*M,N jj[ Jo(N) 
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where g T dg is the Haar measure on O(N) and Zm,n is a normalization constant. The 
j.p.d.f. for complex and quarternionic Wishart ensembles are similar 

P 2 (A) = J-\A(\)\*f[\?- N [ e-^^Vdg, 

Z M,JV jj[ Ju(N) 

for complex Wishart and 

N r 

P 4 ( A ) = _L_ |A(A)| 4 fTA? M - JV)+1 / e- 2 ^^" 1 ^- 1 ))^ 1 ^, 

for the quarternionic Wishart. 

A major difficulty in the asymptotic analysis of the real Wishart ensembles is to find 
a simple expression for the j.p.d.f. in terms of its eigenvalues. For complex Wishart 
ensembles, the integral over U(N) in the expression of the j.p.d.f. can be evaluated using 
the Harish-Chandra [42] Itzykson Zuber [32] formula to obtain a compact determinantal 
formula for the j.p.d.f. in terms of the eigenvalues Xj. This has led to the results in 
[TO] and results on the largest eigenvalue distributions for more general complex Wishart 
ensembles [38], [5?]. In the quarternonic case, although the Harish-Chandra and Itzykson 
Zuber formula does not apply, the integral over Sp(N) in the j.p.d.f. can still be evaluated 
using Zonal polynomial expansion to obtain a determinantal formula for the j.p.d.f. in the 
rank 1 case [73]. In the real case, however, neither of these methods apply and so far there 
has not been any simple formula for the j.p.d.f. that allows the asymptotic analysis in the 
real case. Our first result is a contour integral formula that would allow the asymptotic 
analysis of the rank 1 real Wishart spiked model. 

Theorem 2.1. Assuming N is even. Let the non-trivial eigenvalue in the covariance 
matrix £ be 1 + r and suppose t ^ 0. Then the j.p.d.f. of the eigenvalues in the rank 1 
real Wishart spiked model with covariance matrix X is given by 

/N M-N-l 

where V is a contour that encloses all the points Ai, . . . , Ajv that is oriented in the counter- 
clockwise direction and Zm,n is the normalization constant. The branch cuts of the square 
_ i 

root (t — x) 2 is chosen to be the line arg(t — x) — it. 

We will present two different proofs of this in the paper. The first one is a geometric 
proof which involves choosing a suitable set of coordinates on O(N) and decompose the 
Haar measure into two parts so that the integral in ( 12. ip can be evaluated. This will be 
achieved in Sections 13 . 1 1 and 13 . 2 1 The second proof is an algebraic proof that uses the Zonal 
polynomial expansion to verify the formula in Theorem 12.11 This proof will be given in 
Appendix A where integral formulae of the form ( 12. 2 p for the complex and quarternionic 
Wishart ensembles will also be derived. 
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Remark 2.1. The integral formula derived here is very similar to a more general formula 
in fT3tf . in which the matrix integral over 0(N) is given by 



( w Vd<? oc / — T dS 

o{N) J H j=1 det(S -yjX)2 



e" tr| 



where the integral of S is over y— T times the space of N x N real symmetric matrices and 
yj are the eigenvalues of Y . The measure dS is the flat Lebesgue measure on this space. 
The proof in [13] is also similar to our geometric proof of \2. 2\) presented in the main text. 
The main difference between the derivation in [13] and our derivation is that in obtaining 
h2.2\) . we use the assumption that Y is a rank 1 matrix to decompose the Haar measure to 
obtain the simpler formula \2. 2\) in the rank 1 case. 



Remark 2.2. Shortly after the first part of this paper, which contains h2.2\) appeared in 
the preprint server ArXiv, we learnt that D. Wang has also derived \2. 2\) independently 
and use it for the asymptotic analysis in the super- critical regime [71 



From the expression of the j.p.d.f., we see that the largest eigenvalue distribution is 
given by 



P(A ma ,. <z)= I ... f P{X)dX x . . . d\ N , 

J \ 1 <...<\ N <z J 

= Zj? N [ e^j f ... f \A(X)\T]w(X j )dX 1 ...dX N dt 



(2.3) 



M _ M-N-1 



where w(x) is 

w(x) = e~^~ x x"'~2^ (t - x)~ f (2.4) 

and T is a close contour that encloses the interval [0, z]. 

We can analyze the integrand as in [53], [68], [67] and [69J. By an identity of Brujin 
we can express the multiple integral as a Pfaffian. 



N 



A(A)in^(Ai)dA 1 ...dA JV 
\ 1 <...<x N <z J j=1 (2.5) 

= p f «(! - X[z,oo))rj(x), (1 - X[z,oo))rk{y))^ ■ 

where Tj(x) is an arbitrary sequence of degree j monic polynomials and {f,g) 1 is the skew 
product 

poo />oo 

(f,g)i= / e(x -y)f(x)g(y)w(x)w(y)dxdy. (2.6) 
Jo Jo 

where e(x) = |sgn(x). In defining the skew product, the contour of integration will be 
defined such that if t is on (0,oo), then the interval (0, oo) will be deformed into the 
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upper or lower half plane near t such that the integral is well defined. As T will only 
intersect (0, oo) at a point xq > z, the left hand side of (12.51) . and hence the Pfaffian, 
will not be affected by such deformations. Then by following the method in [68], [67] 
and [6H], we can write the Pffafian as the square root of a Fredholm determinant. Let 
D be the moment matrix with entries {Tj,r k ) 1 and det2 be the regularized 2-determinant 
det 2 (J + A) = det((J + A)e~ A )e tr{ - All+A22) for the 2 x 2 matrix kernel A with entries A ij} 
then we have 

Pf «(1 - X[z,oo))(x) rj (x), (1 - Xiz^Mniy))^ = VdetD^Vdeta {I- X K X ), 
where % = X\z,oo) an d K is the operator whose kernel is given by 

K(xv)=( Sl{X,y) -i Sl ^ y) ) (2 7) 

and Si(x,y) and ISi(x,y) are the kernels 

N-l 

Si(x,y) = -^2 r j (x)w(x)fj ljk e(r k w)(y), 
ISi(x,y) = t(rjw)(x)/j, jk e(r k w)(y) 

j,k=0 

and fij k is the inverse of the matrix D. As shown in [73J , the kernel can now be expressed in 
terms of the Christoffel Darboux kernel of some suitable orthogonal polynomials, together 
with a correction term which gives rise to a finite rank perturbation to the Christoffel 
Darboux kernel. In this paper, we introduce a new proof of this using skew orthogonal 
polynomials and their representations as multi-orthogonal polynomials. By using ideas 
from [2J to write skew orthogonal polynomials in terms of orthogonal polynomials, we can 
express the skew orthogonal polynomials with respect to the weight w(x) in terms of a 
sum of Laguerre polynomials. Let Tr kj i be the monic skew orthogonal polynomial of degree 
k with respect to the weight w(x). 

<vr 2fc+ i,i, y j ) 1 = (tt 2M , yi) l = 0, j = 0, . . . , 2k - 1. (2.9) 

Then we can write these down in terms of Laguerre polynomials. 

Proposition 2.1. Let L k be the degree k monic Laguerre polynomial with respect to the 
weight w (x) 

/>oo 

L k (x)Lj(x)w (x)dx = 5j k hjfl, w (x) = x M ~ N e~ Mx . 

If (Z/2fc-i, £2*1-2)1 7^ 0> then the skew orthogonal polynomials iT2k,i o>nd'K2k+i,i both exist and 
Tr 2k 1 is unique while vr 2 fc + i 5 i is unique up to an addition of a multiple ofir^ki- Moreover, 
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we have (L 2 k, L 2 k-i)i = and the skew orthogonal polynomials are given by 



^2k,l — J^2k ~ TT 7 r--^2fe-l, 

W2t-1, -^2fc-2/i 

T (L 2 k+l, L 2 k-2)i T \L 2 k+li L>2k-l)i T 

7I"2fc+l,l — J->2k+l — TT j r--^2fc-l + 77 j ^—J-'2k-2 + C' K 2k,li 

W2fc-1, -^2fc-2/i W2fc-1, -^2fc-2/i 

where c is an arbitrary constant. 

Next, by representing skew orthogonal polynomials as multi-orthogonal polynomials 
and write them in terms of the solution of a Riemann-Hilbert problem as in [55] . we can 
apply the results of [28] and [9] to express the kernel Si(x,y) as a finite rank perturbation 
of the Christoffel Darboux kernel of the Laguerre polynomials. 

Theorem 2.2. Suppose (Ljv-ij L N _ 2 ) 1 (L N _ 3 , L N _ 4 ) 1 ^ 0. Let S\(x,y) be defined by K2. 8\) 
and choose the sequence of monic polynomials Tj(x) such that rj(x) are arbitrary degree j 
monic polynomials that are independent on t and rj(x) = 7Tj^(x) for j = N — 2, N — 1. 
Then we have 



Si(x,y) - K 2 (x,y) = 

e (n N+hl w ir Ntl w) (y) [ M Mt-ytf+N) ] ( 2nTi(x) ) w(yX 



-2H^To \ (L N . 2 {x)\ , , (2.10) 



2/lJV-2,0 2/ljv-l,0 / 

where K 2 (x,y) is the kernel of the Laguerre polynomials 



i 



Note that the correction term on the right hand side of (12.101) is the kernel of a finite rank 
operator. Its asymptotics can be computed using the known asymptotics of the Laguerre 
polynomials and the method in [S3], [S] and [32] • The actual asymptotic analysis of 
this correction term, however, is particularly tedious as one would need to compute the 
asymptotics of the skew orthogonal polynomials up to the third leading order term due 
to cancelations. To compute the contribution from the determinant detO, we derive the 
following expression for the logarithmic derivative of det D. 

Proposition 2.2. Let D be the moment matrix with entries {Tj,rk) v where the sequence of 
monic polynomials rj(x) is chosen such that rj(x) are arbitrary degree j monic polynomials 
that are independent on t and rj(x) = TZj \{x) for j = N — 2, N — 1. Then the logarithmic 
derivative of det D with respect to t is given by 

|-logdetD = - I Sli - X ' X h x. (2.12) 
at L , t - x 
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This then allows us to express the largest eigenvalue distribution F(X max < z) as an 
integral of determinant. 

Theorem 2.3. The largest eigenvalue distribution of the rank 1 real Wishart ensemble can 
be written in the following integral form. 

P(A ma , < z) = C [ exp ( -Ml* ~\f I ^^dxd S \ v/det 2 (/ - X K X )dt. 

JT VA i + r J z Jc JR + S — X J 

(2.13) 

for some constant Cq and K is the operator with kernel given by Ii2.8\) . The integration 
contour V is a close contour that encloses the interval [0, z] in the anti- clockwise direction. 

These results so far are valid for all N and M and are exact. As N, M — > oo, the 
integral expression (12. 13p can be used for the asymptotic analysis to obtain the largest 
eigenvalue distribution when the phase transition occur. In this paper, we will demonstrate 
the asymptotic analysis for the case where M/N — > r y 2 = 1, but the analysis can also be 
applied to the case of any finite 7. 

In the asymptotic limit, the t integral in (I2.13P can be computed using steepest descent 
analysis. In fact, we shall see that when t ^ b + , the integrand in (12.131) is of order 



/ Mrt 1 /"* f SAx^x) , , \ n — „ 

exp — -nil dxds Vdets (7 - X K X ) ~ Ce^ 

\2(l + r) 2J J R s-x J 



'hSm + P(s) log(t-s)ds 



where p(s) is given in (II. ip . The saddle point t sa ddie of this equation does not belong to 
the bulk region [&_, b + ] unless r is at the critical value r c = 7^. When r < r c , the steepest 
descent contour T can be deformed such that it does not intersect [&_,&+.]. In this case, t 
will always be of finite distance from [&_, b + ] and the factor (t — x) _ 2 in w(x) will have no 
effect on the asymptotics of the kernel Si at the edge point b + . In this case, the kernel at 
the edge point will be given by the Airy kernel. 

2 2 
where £1 = (x — 6+) 7M \ and £2 = ("V — ^+) 7M ^4 are finite. 

(1+7)3 (1+7)3 

In the critical case, however, the main contribution to the contour integral (I2.13P comes 
from a small neighborhood of t = b + and the factor (t — x)~z in w(x) will now significantly 
alter the behavior of the kernel Si and change it from the Airy kernel to into a more 
complicated kernel. This gives rise to the phase transition and a new distribution function. 
Our next result is the representation of this new distribution function in terms of the 
Hastings-McLeod solution of the Painleve II equation. First let us define some functions 
that will appear in our formula. 
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Let ip(u,T) = (T — u)2 and Hj(u,T) = Ai^\u)ip 1 (u,T) and let «Sio, <Sn be 



Hi(u, s)du + 



oo 
oo 



Ht / H 2 dv 



S 10 (T) = 

Sn(T) = 2 

+ 2 J E x I H 2 dvdu 
Define the function S(T) to be 



(-u)s 



2^3 

-d s --T,. 



dw 



2tt(T - u) 

#idw / H 2 du. 



pT poo 

S(T)=S 10 (T)+ 5u(s)dr-l/2 / (5u j+ (s)-5i 1 ,_(s))ds, (2.14) 
Jo «/ 

where the contour of integration in the first term remains in the upper half plane and <Sii,± 
are the boundary values of Sn as it approaches the real axis in the upper/lower half plane. 
Now let U be the matrix 



U 



( 






-7/^0 \ 

^Vo -V>~Vo 



J 



1 

2i/> 2 



and define /ij to be the vector hj = yO, 0, 0, ^-J for j = 0, 1, 2 and /ij = for j = 3 and 
j = 4, where O is the Hastings-McLeod solution of Painleve II (11.31) and 



*(C) 



n + adlr 



0) 



02= [C + J MZ)MZ)*Z)<Po-<r<h- 

Let u,- be the vector that satisfies the linear system of ODEs with the following boundary 
condition 



dVj_ 



U(C)vj + hj, vj ~ 0, 0, 0, / Hjdu , C ->• +00, j = 0, 1, 2 



(2.15) 



t? 3 ~(0, 0,-1,0), C^-oo, tf 4 ~(0, 0,-1,0), C^+oo 
Then the largest eigenvalue distribution at the phase transition is given by 
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Theorem 2.4. Suppose N is even. Let w = ("j) 3 (1 — r) G (—00,00) and let ( = 

(z — 4) (iV/4) 1 , then as N, M — » 00 snc/i £/ia£ M/iV — >• 7 2 = 1, the largest eigenvalue 
distribution at the phase transition is given by 



Km P f(A ma:E - 4) f ^ 3 < = CyTjy 2 (C)Im 




x (det (5 jfc - (a.j,/3k)) 



i<i,fc<3 




/or some constant C , where S + zs a contour in the upper half plane that does not contain 
any zero of f_ H$du and approaches 00 in the sector 11/ 3 < argT < n. It intersects R at 
£/je point £. The entries in the 3x3 matrix are given by 

v 34 - v 33 + 1 ( a . v 32 ( n ^ ^ w («33, 
(ai,pi) = , (ai,/J 2 ) = — , (ai,Pi 



2 ' v ±,r ^ 2a 
T ~ 1 

(a 2 , p\) = —go - Bigi - <?2 - -B 2 ^_, (a 3 , ft) = -gi + i3ig 

T 

(a 2 , = — - &1U-M-1 - M_ j2 j-2 - B 2 V- s j-2, 

(a 3 , = ^ M -,ii-2 + #1^-2,2, j = 2, 3, 

where W(f,g) is the Wronskian W(f,g) = fg' — gf and Vjk are the components of the 
vectors Vj in 112.15}) . The functions B\, B 2 , qj, R-, P_j and u~jk are given by 

qj = v j4 - v j3 - J Hjdu, u-jk = , J = 1, 2, 



7£_ = W44 - ^43 + 1, V-fi = V 42 , V-A = 



a 

00 



* X 2/.>,(«)d U ' 

B 2 = -^-- H du + ^- H 1 du + - H 2 du. 

* J —oo " J —oo " J —oo 

The distribution in Theorem 12.41 is expressed in terms of integrals of the Airy function, 
together with the functions Vjk, which are solutions of linear ODEs with known boundary 
conditions. The coefficients of the ODEs satisfied by the Vjk are given in terms of the 
Hastings-McLeod solutions and its derivatives and are therefore known functions. 



Remark 2.3. The distribution in Theorem \2.J\ can also be expressed in terms of an in- 
tegral involving the solution of a Riemann-Hilbert problem, (see Appendix B) With the 
recent advancement in the numerical computation of Riemann-Hilbert problems [56], this 
representation may be useful for the numerical computation of the distribution in Theorem 
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2.1 Other results 



In obtaining the main result in Theorem I2.4[ we have obtained some other results which 
may also be of interest to mathematicians and physicists working in random matrix theory. 

2.1.1 Random matrix with external source 

One applications of the orthogonal polynomial approach developed in this paper is in the 
studies of random matrix with a rank 1 external source. Random matrix with external 
source are random matrix models on the space of real symmetric, Hermitian or Hermitian 
self-dual N x N matrices with the following probability measure 

z 

for some real-valued function V(x) such that e~ v ^ decays fast enough as x — > ±oo and a 
real symmetric, Hermitian or Hermitian self-dual matrix A. The function V(x) is usually 
called the potential while the matrix A is known as the external source. The real Wishart 
ensemble can be thought of as a special case when V(x) = f — "jl" 1 logx. Random 
matrices with external source are first studied by Brezin and Hikami [25], [26] and P. 
Zinn- Justin [71], [75] as a model of systems with both random and deterministic parts. 
For Hermitian random matrices, the external source model can be studied using multi- 
orthogonal polynomial and Riemann-Hilbert techniques JTS] and there are many recent 
advancements in the asymptotic analysis of Hermitian external source models with non- 
Gaussian potential V{x) [32], [H], [32], 02], [20], [21], [22]. In [72], the contour integral 
formula (12. 2p was derived independently and was used to study real symmetric random 
matrix with a rank 1 external source and a large class of potential V(x). By using the 
linear statistics results of Johansson [UJ and a representation for the largest eigenvalue 
distribution that is different to (12.131) . the largest eigenvalue distribution was obtained in 
the super-critical regime. The orthogonal polynomial approach developed in this paper can 
be used to extend the results in [72J to the critical regime. In fact, for real random matrix 
with a rank 1 external source, the expression (12.131) for the largest eigenvalue distribution 
remains valid, although the kernels Si and K will have to be modified according to the 
potential. For a polynomial potential V(x), the analogue of (I2.1(jp . which expresses the 
kernel S\ in terms of orthogonal polynomials is well-known [73] . By using the asymptotics of 
orthogonal polynomials found in [36] and [37], such representation can be used to compute 
the asymptotics of the kernel Si, which can then be used in ( 12 . 1 3[) to obtain the largest 
eigenvalue distribution. 

2.1.2 Orthogonal ensembles 

A large part of this paper involves the analysis of an orthogonal ensemble with weight (12 .4p 
and in doing so, we have obtained some new results for orthogonal ensembles. 
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The first of these results is the logarithmic derivative of the partition function Z\ = 
det D in Proposition 12.21 While the proposition is stated in the terms of the derivative 
of the parameter t in w(x), the proof can easily be generalized to obtain the logarithmic 
derivative of det D for a general orthogonal ensemble with respect to any parameter. 



For polynomial potential, the asymptotics of Si(x, x) can be obtained through the asymp- 
totics of orthogonal polynomials found in [36J, [3?]. This could then be used to compute 
the asymptotics of the partition function Z\ = detD. (For unitary ensembles, this was 
done in [T7| for the quartic potential V(x) — \ + ^- using a different differential identity) 
Asymptotic analysis of the partition function is of importance in extending the universal- 
ity results of orthogonal and symplectic ensembles to general weights w. At the moment, 
universality in the orthogonal and symplectic ensembles are proven for a large class of po- 
tentials V(x) [31], [33], [63], [60], [61]. However, except for the quartic case V(x) = \ + %-, 
all the available results are restricted to the case where the limiting eigenvalue distribution 
is supported on a single interval. The main obstacle in extending these results to more 
general potential is the computation of 



in the limit N — > oo, where Z n ^{w) is the partition functions of the ensembles (See remark 
2.4 of [63] and remark 1.5 of [3~T] ) 



In order to proof the universality in orthogonal and symplectic ensembles using the method 
in [21], [22] and [63J, one needs to show that limjv^oo det Q ^ 0. While the leading order 
asymptotics of these partition functions can be found using the estimates in [47], their 
combined contributions to det Q cancel in the leading order and hence higher order terms 
in the asymptotics of Z n ^{w) are needed to show that liniAr_ 5 . 00 det Q ^ 0. At the moment, 
asymptotics for the sub-leading order terms of partition functions are only available for 
(5 = 2. By using a differential identity for log Z N<2 , the authors in [T7J computed the sub- 
leading order terms in Zn,2 for the potential V(x) = ^ + ^- as iV — > oo. A combination 
of the method in [T7j and the differential identity (12.161) may provide a way to compute 
the sub-leading order terms in the partition function Z^,i and help extend the universality 
results in orthogonal and symplectic ensembles to more general potentials. 

Another interesting observation is Corollary 14. 11 in which we showed that (L^, i^fc— 1)1 — 
whenever k is even. This turns out to be a very useful identity in the analysis of the 

2 

phase transition when the double scaling limit t — 4 = T(4/iV) 3 for the ensemble with 




(2.16) 
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weight ( 12 .4p has to be considered. When analyzing this double scaling limit, the identity 
(Lk, Lk-i) 1 = leads to cancelation in the leading order terms of the kernel Si. This 

2 4 

enables us to show that S\ is of order N3 instead of A^3, which is essential for the scaled 
limit of the determinant det2 (/ — X^-X) to exist. We believe this type of identity will also 
be useful in the analysis of other double scaling limits in orthogonal ensembles. 

The paper is organized as follows. In Section [3] we will prove the contour integral 
formula (12.21) and in Section H] the Christoffel-Darboux formula (12.101) for the kernel Si will 
be derived. In Section |5l we will prove the differential identity for the moment matrix D 
in Proposition 12.21 The results in these sections are all exact and apply to all N and M. 

We will start the asymptotic analysis in Section |6] in which the asymptotics for the kernel 
Si, det D and det2 (I — X^-X) wm De obtained. Finally, we will express the asymptotics of 
the determinant det2 (/ — X^X) i n terms of the Painleve transcendents in Section [7J 

Throughout the paper, we shall assume that N is even and that M — N > 0. 



3 Contour integral formula for the j.p.d.f. 

In this section we will prove the integral formula for the j.p.d.f. in Theorem 12.11 



3.1 Haar measure on SO{N) 

In this section, we will find a convenient set of coordinate on O(N) to evaluate the integral 



that appears in the expression of the j.p.d.f. (12.11) . As both E -1 and Y are symmetric 
matrices, they can be diagonalized by matrices in O(N). We can therefore replace both 
and Y by the diagonal matrices S^" 1 and A<j. 



The group O(N) has two connected components, SO(N) and GL(iV) that consists of 
orthogonal matrices that have determinant 1 and —1 respectively. Let T be the matrix 
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then the left multiplication by T defines an diffeomorphism from O-(N) to SO(N). In 
particular, we can write the integral over O(N) in (12. ip as 



I(£,A)= / e-f^" 1 ^ 1 )/^, 

JsO(N) JO-(N) 



(N) JO.(N) 



SO(iV) JSO{N) 

= [ e-f t ^ A ^g T dg+ [ e-^^ A ^g T dg, 
where is the diagonal matrix with the first two entries of swapped. 



~_, ,. / 1 1 1 \ 
L . = diag , . . . , . 



Note that g T dg is also the Haar measure on SO(N). 

As we are considering the rank 1 spiked model, we let T\ — . . . — 7>/_i = and r N = r. 
Therefore = and we have 

J(E, A) = 2 / e-f^W- 1 )^ ( 31 ) 

JSO(N) 

Let <7y be the entries of (7 G SO(N). Then the integral J can be written as 



/(£, A) = 2 I e-f tr ^ l9A ^g T dg, 

JsO(N) 



>(JV) 

e -f tr ((s^ 1 -/iv)9A d 3- 1 ) e -ftr( 9 A d9 - 1 )^T^ ) 



SO(N) 



N . 

j=l Jso(N) 



We will now find an expression of the Haar measure and use it to compute the integral 
J(E,A). 

First let us define a set of coordinates on SO(N) that is convenient for our purpose. 
We will then express the Haar measure on SO(N) in terms of these coordinates. 
An element g G SO(n) can be written in the following form 



9 = (9i, ■ ■ ■ , 9n) , \9i\ = 1, 9i ■ 9j = S. 



1 j ■ 



This represents SO(N) as the set of positively oriented orthonormal frames in M. N whose 
coordinate axis are given by the vectors <fj. As the vector gpf is a unit vector, we can write 



16 



its components as 

3-1 

g 1N = cos(j) 1 , g jN = JJsinfocos^, j = 2, . . . , n - 1, 

fc=1 (3.2) 

N-l v y 

9nn = Y[ sin 

fe=i 

The remaining vectors gx, ■ ■ ■ ,Qn-i form an orthonormal frame with positive orientation 
in a copy of M^ -1 that is orthogonal to g^. Therefore the set of vectors g%, . . . ,g*N-i 
can be identified with SO(N — 1). To be precise, let u be a unit vector in M. N and let 
G(u) e SO(N) be a matrix that maps w to the vector (0, . . . , 0, 1) T . As G is orthogonal, 
we have 

rp 

G{gN)gj = (vji, 0) , j < N (3.3) 

In particular, the matrix V with entries for 1 < z, j < iV — 1 is in SO(N — 1). The 
following then gives a set of coordinates on SO(N). 

9 = (9N,V). (3.4) 

In the above equation, pjv is identified with the coordinates <pj in (13.21) . while the matrix 
V is identified with the coordinates in SO(N — 1) that correspond to V. In terms of these 
coordinates, the left action of an element S G SO(N) on g is given by the following. 

Sg = {Sgi, Sg N -!, Sg N ) T 

= (SG(g N )~ 1 Vi, SG^gN^VN-!, Sg N ) T 

Then as in (I3.3p . we have 

G{Sg N )SG{g N )~ 1; Vj = . . . , %,v-i, 0) T . 

The matrix V with entries are again in SO(N-l), therefore the matrix G(SgN)SG(gi^) 
is of the form 

GiSg^SGifay 1 = g Q (3.5) 

From the fact that GlSg^SGlg^)" 1 is an orthogonal matrix, it is easy to check that 
s = and = ±1. To determine sjy, let us consider the action of G^Sg^SG^giy)^ 1 on 
(0, 0, ... , 1) T . We have 

G^SG^vr^O, 0, . . . , If = G(Sg N )Sg N = (0,0,..., if 

Therefore = 1 and Sjv-i is m SO(N — 1). The action of 5 on g is therefore given by 

Sg= (s$ N ,S N -iVy (3.6) 
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We will now write the Haar measure on SO(N) in terms of the coordinates (13. 4p . These 
coordinates give a local diffeomorphism between SO(N) and S 1 ^ 1 x SO(N — 1) as G 
S^ -1 and V G SO(N — 1). Let c?X be a measure on S N_1 that is invariant under the action 
of SO(N) and V T dV be the Haar measure on SO(N — 1), then the following measure 

dH = dX A V T dV, 

is invariant under the left action of SO(N). Let 5 G SO(N), then its action on the 
point (<fjv, V") is given by ()3.6|) . where Sjv-i depends only on the coordinates (pi, ... , 4>n-i- 
Therefore under the action of S, the measure dH becomes 

dH dX A y T ^_ 1 5 7V -i^ = c/X A y T dV, (3.7) 

as <iX is invariant under the action of S. Therefore if we can find a measure on S^ -1 that is 
invariant under the action of SO(N), then dX A V T dV will give us a left invariant measure 
on SO(N). Since the left invariant measure on a compact group is also right invariant, this 
will give us the Haar measure on SO(N). As the metric on S 1 ^' 1 is invariant under the 
action of SO(N), it is clear that the volume form on S' 7V_1 is invariant under the action of 
SO(N). Let dX be the volume form on S^ -1 , then from (13. 7J , we see that the measure 
dX A V T dV is invariant under the action of SO(N). 

Proposition 3.1. Let dX be the volume form on S 1 ^^ 1 given by 

dX = sm N ~ 2 ((/>{) sin^- 1 ^) . . . wn((f> N - 2 ) /\f~i dfa 



in terms of the coordinates <f>i, . . . , 0jv-i i n A3.2\) and j3.4ty , then the Haar measure on 
SO(N) is equal to a constant multiple of 

dH = dX A V T dV, 

where V T dV is the Haar measure on SO(N — 1) in terms of the coordinates 
We can now compute the integral /(£, A). 

3.2 Integral formula 

By using the expression of the Haar measure derived in the last section, we can now write 
the integral /(£, A) as 

N . 

/(S, A) = 2TTe~~ J / e^+^^= lX ^ N g T dg, 
j=l Jso(N) 

N . . 

= 2T[e-% Xi / V T dV / e^)^= lX ^ N dX, 

■ =1 JSO(N~l) JS"- 1 

N „ 

e --^ / e w+^^ x ^dX, 

: , JS N ' 1 



18 



for some constant C, where the N — 1 sphere S 1 ^ -1 in the above formula is defined by 
J2f=i 9jN = 1 an d dX is the volume form on it. If we let gjjy = Xj, then the above can be 
written as 

N n / N \ 

J(E, A) = 2C ]7e ^ / e wih)^U x ^s [^x] - l\ dx x . ..dx N . (3.8) 
i=i JrN W J 

This can be seen most easily by the use of polar coordinates in R , which are given by 

o-i 

X\ — r cos 0i, Xj = r J^J sin0 fc cos0j, j = 2,...,N—l, 

k=l 

N-l 

xn = r n sin ^' 

k=l 

Then the volume form in ~R N is given by 

dx\ . . . dxN = r N ~ l sin^ -2 0i . . . sin (f)N^ 2 drd(f)i . . . d<pN-i 
Therefore in terms of polar coordinates, we have 



"N 



noo ( N \ 



7T p2tv 



o </o 

tAI sr^N \ 2 „ 

x e^+^ J ^ sin^" 2 0i . . . sin 07v_ 2 



tai sr N \ ^2 
e2(1+T) |A '-r/.Y. 



S N-1 



To compute the integral J(S, A), we use a method in the studies of random pure quantum 
systems (see, e.g. [52]). The idea is to consider the Laplace transform of the function 
I(E,A,t) defined by 

n M n ( N \ 

/(£, A, t) = 2C TT e"^ / e^T E ^ AjX ?5 V a; 2 - t \ dx x . . . dx N , 

f=i J ^ N \U J 

then /(£, A, 1) = /(£, A). The Laplace transform of /(£, A, t) in the variable t is given by 



/•OO J v /* 

/ e- st I(Z,A,t)dt = 2CTTe-^ / 
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Then, provided Re(s) > max.,- (Xj), the integral can be computed explicitly to obtain 



Taking the inverse Laplace transform, we obtain an integral expression for J(E,A). 



where T is a contour that encloses all the points 2 (x+ T ) ' ' ' ' 2(i+rj ^ ^at * s or i en t;ed in 
the counter-clockwise direction. Rescaling the variable s to s = Mt, we obtain 



This then give us an integral expression for the j.p.d.f. in Theorem 12.11 

For the purpose of computing the largest eigenvalue distribution F (\ max < z), we can 
assume that the eigenvalues are all smaller than or equal to a constant z. 

4 Skew orthogonal polynomials and the kernel Si 

In this section we will prove the representation of the kernel Si in Theorem 12.21 We will do 
so by using the multi-orthogonal polynomial representation of skew orthogonal polynomials 
in [70] and then apply the Christoffel-Darboux formula for multi-orthogonal polynomials 
in [28] to write the kernel Si as a finite sum of multi-orthogonal polynomials. We then 
simplify this sum further by using a result in [9j. This gives a new proof to a well-known 
result of Widom [73]. 

4.1 Skew orthogonal polynomials 

As explain in the introduction, we need to find the skew orthogonal polynomials with the 
weight (12. 4p . Let us consider the skew orthogonal polynomials with respect to the weight 
w in (12. 4p . We shall use the ideas in [2] to express the skew orthogonal polynomials in 
terms of a linear combinations of Laguerre polynomials. 
Let Qj{x) to be the degree j + 2 polynomial 






Qj(x) 



d 

dx 



(x J+ (t — x)w{x)j w 1 (x), j > 0. 



Then as we assume M — N > 0, it is easy to see that 




(4.1) 
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for any f(x) such that J °° f(x)w(x)dx is finite, where the product () 2 is defined by 

POO 

(f(x)g(x)) 2 = / f(x)g(x)w (x)dx, w (x) = x M ~ N e~ Mx . (4.2) 



Note that w (x) is not the square of w(x). The fact that w (x) is the weight for the 
Laguerre polynomials allows us to express the skew orthogonal polynomials for the weight 
( 12. 4 p in terms of Laguerre polynomials. 

In particular, this implies that the conditions ( 12. 9 p is equivalent to the following con- 
ditions 

<vr 2M ,l/ i > 1 = 0, J = 0,1, 
<7r2M^) 2 = °> j = 0,...,2fc-3. 

and the exactly same conditions for 7T2fc+i,i(x). In particular, the second condition implies 
the skew orthogonal polynomials can be written as 

vr 2 fc,i(x) = L 2k (x) + 72fc,i^2fc-iO) + j2k,2L 2k - 2 (x), 

7T2fe+l,l(^) = L 2 k+l{x) + 72fe + l,0^2fe(^) + 72fe+l,l-^2fc-l(^) + 72fe+l,2-^2fe-2 {%) , 

where Lj(x) are the degree j monic Laguerre polynomials that are orthogonal with respect 
to the weight wq(x). 

Ln{x) = [ l > 6 X ±_ r e -M* n+M-Nj 

= x"- (M ~ iV + " ) V- 1 + 0(x"- 2 ). 

M v ' 

The constants j^j are to be determined from the first condition in (14.31) . We will now show 
that if (L 2 k-i, £2*1-2)1 7^ 0> then the skew orthogonal polynomials TT 2 k,i an d ^2^+1,1 exist 
and that ~K 2 k,i is unique. First let us show that the first condition in H4. 3j) is equivalent to 

(v"2fc,i, L 2 k-j) 1 = 0, j = 1, 2. 

To do this, we will first define a map qn from the span of L 2 k-i and L 2 k- 2 to the span of 
y and 1. 

Let P(x) be a polynomial of degree m. Then we can write the polynomial P(x) as 

P(x) = — (q(x)x(t — x)w(x)) w^ 1 (x) + R(x) (4.5) 
ax 

where q(x) is a polynomial of degree m — 2 and i?(x) is a polynomial of degree less than 
or equal to 1. By writing down the system of linear equations satisfied by the coefficients 
of q(x) and R(x), we see that the polynomials q(x) and R(x) are uniquely defined for any 
given P(x). In particular, the map / : P(x) 1— > R(x) is a well-defined linear map from the 
space of polynomial to the space of polynomials of degrees less than or equal to 1. Let Qk 
be the following restriction of this map. 
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Definition 4.1. For any polynomial P(x), let f be the map that maps P(x) to R(x) 
in ( |^.5| ). Then the map g k is the restriction of f to the linear subspace spanned by the 
orthogonal polynomials L k , L k _\. 

We then have the following. 

Lemma 4.1. // (L k , L k -\) l 7^ ; then the map g k is invertible. 

Proof. Suppose there is exists non-zero constants a\ and 02 such that 

a\L k + a 2 L k _i = -r(q(x)x(t — x)w)w^ 1 
ax 

for some polynomial q(x) of degree k — 2, then by taking the skew product () 1 of this 
polynomial with L k , we obtain 

a 2 (L k -x, L k ) 1 = (a x L k + a 2 L k - X , L k ) x = (q(x), L k ) 2 = 0. 

As q(x) is of degree k — 2. Since (L k _i, L k ) 1 7^ 0, this shows that a 2 = 0. By taking 
the skew product with L k -i, we conclude that a\ = and hence the map g k has a trivial 
kernel. □ 

In particular, we have the following. 

Corollary 4.1. If k is even, then (L k , L k ~i) 1 = 0. 

Proof. Let q(x) be a polynomial of degree k — 2 that satisfies the following conditions 

— (q{x)wi{x))x^Wi{x)dx = 0, j = 0,...,k — 2, (4.6) 

where Wi(x) = x 2~^(£— x)^e~^ . A non trivial polynomial q(x) of degree k— 2 that satis- 
fies these conditions exists if and only if the moment matrix with entries J R -^(x l w^{x))x^Wi{ 
has a vanishing determinant. For even k, the moment matrix is of odd dimension and anti- 
symmetric and hence its determinant is always zero. 

Assuming k is even and let q(x) be a polynomial that satisfies (14. 6p . By taking the 
inner product () 2 with x- 7 , we see that there exists non-zero constants a\ and a 2 such that 

aiL k + a 2 L k _i = —{q{x)x{t — x)w)w^ 1 , 

Therefore by Lemma [4. 11 we see that if k is even, we will have (L k , L k _i) l = 0. □ 

Lemma I4TT1 shows that if (Li, Lj_i) 1 7^ 0, then there exist two independent polynomials 
Ro(y) and Ri(y) in the span of y and 1 such that Rj(y) = g^L^j). Then we have 

Rj{y) = —j-(qj(y)v(t - y)w)w~ l + L^^y), j = o, 1. 
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In particular, the skew product of Rj{y) with I < 2 is given by 

(Li_i(x), R j (y)} 1 = - (Li-i, qj} 2 + (Li-i, L i _ j ) 1 . 

As qj is a polynomial of degree less than or equal to i — 2 and I < 2, the first term on the 
right hand side is zero. Therefore we have 

(Li-i(x), Rj(y)) l — (Li_i, L i ^j) l , / < 2, j = 0, 1. (4.7) 

We can now show that the skew orthogonal polynomials ir 2 k,i and 7T2ft+i,i exist if (L 2 k-i, £2*1-2) 
0. 

Proposition 4.1. If (L2k-i, -^2fc-2)x 7^ 0, i/ien £/ie s£;eu> orthogonal polynomials 7^,1 and 
7T2fe+i,i e:risi and 7T2fc,i unique while TT 2 k+i.i is unique up to an addition of a multiple 
of 7T2fe,i- Moreover, we have (L^k, ^2fc-i)i — ^ e sfceu> orthogonal polynomials are 
given by 

T (Lzk-i -^2fc-2/i r 

7T2fc,l — -^2fc — TT 7 T--^2fe-l, 

\^2fc-l,^2fc-2/i ^ 

r (-^2/c+l, L 2 k-2)i T (L2k+1, L 2 k-l) 1 r 

7T2fe+l,l — -t>2fc+l — TT 7 T--^2fc-l + 77 7 T--^2fc-2 + C7T2fc,l, 

W2fc-1, -^2fc-2/i W2A;-1, -^2fc-2/i 

/or k > 2, where c is an arbitrary constant. 

Proof. Let 7r 2fc l and 7^+1,1 be polynomials defined by 

7r2fc,l(x) = L 2 fc(x) + 72fc,l-^2fc-l(x) + 72fc,2-^2fc-2(x), 

7T2jfe+i,i(a;) = L 2 fe+i(x) + 72fc+i,i-^2fc-i(a;) +72^+1,2^-2^), 

for some constants 7^. If we can show that (7r 2fc _/ jl , y 3 '^ = for j = 0, 1 and Z = —1,0, 
then 7r 2 fc— z,i will be the skew orthogonal polynomial. Let i?o an d Ri be the images of 
L 2 k-i and L 2 k- 2 under the map Q 2 k-i- Then by the assumption in the Proposition, they 
are independent in the span of y and 1. Therefore the conditions (iT2k-i,i,y : ') 1 = are 
equivalent to (ir 2 k-i t i, Rj(y)) l = 0. By taking i = 2k — 1 in (J4TTJ) , we see that this is 
equivalent to (7T2k-i,i, L>2k-i-j) 1 = 0. This implies 

\^2k,li L 2 k-l) l = \L2k1 -^2fc-l)i + 72fc,2 \L2k-2, -^2fc-l)i = 0, 
\7T2fc,l) L 2 k-2)i = \L2k7 -^2fc-2/i + 72fc,l \L2k-l, -^2fc-2/i = 0. 



Hence we have 



(L 2 k, L 2 k- 2 ) 1 _ \L2k1 L*2k-i)i 
72fc,i — -77 7 r-, 72fc,2 — 77 7 r~, 



(-^2fc-l, ^2k-2]\ ' (L 2 k~l, L 2 k- 2 ) 1 

which exist and are unique as (L 2 fc_i, L 2 / c _ 2 ) 1 7^ 0. This determines 7^,1 uniquely. By 
Corollary 14.11 we have (L 2 fc, L 2 fe_i) 1 = and hence 7 2 fc )2 = 0. Similarly, the coefficients for 
7r 2 fc+i,i are 

(-^2/e+l, L 2 k-2)i _ (L 2 k+1, L 2 k-\) 1 

72fc+l,l — — 77 7 T", 72fc+l,2 — 77 7 r~- 

\-^2fc-l, -^2fc-2/i \^2fc-l,^2t-2/i 
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Again, these coefficients exist and are unique. However, as (7^,1, ^,1)1 = an< ^ (^fe.i, = 
for j = 0, . . . , 2k — 1, adding any multiple of 7^,1 to 7^+1,1 will not change the orthog- 
onality conditions (7^+1,1, y^) 1 = that is satisfied by 7^+1,1 and hence 7r 2 fc+i ) i is only 
determined up to the addition of a multiple of r K2k ) \- D 



4.2 The Christoffel Darboux formula for the kernel 



In [52], skew orthogonal polynomials were interpreted as multi-orthogonal polynomials and 
represented as the solution of a Riemann-Hilbert problem. This representation allows us 
to use the results in [28] to derive a Christoffel-Darboux formula for the kernel (12. 8p in 
terms of the Riemann-Hilbert problem. 

Let us recall the definitions of multi-orthogonal polynomials. First let the weights Wq, 
Wi and w 2 be 

w (x)=x M - N e- Mx , Wl (x)=w(x) [ e(x-y)L N _ l _ 2 (y)w(y)dy, 1 = 1,2. 



Note that the weights wi(x) are defined with the polynomials L N _ 3 and L N _ 4 instead of 
Ln~i and This is because the construction below involves the polynomial 7Tjv_2,i as 

well as 7Tjv,i- By taking i = N — 3 in (14.71) . we see that the orthogonality conditions for 
7Ttv,i is also equivalent to 

(71^,1,^)2 = 0, j = 0, ...,iV-3, (71^,1,^-^)2 = 0, j = 3,4, 

provided {L N _ 3 , L N _ 4 ) 1 is also non-zero. 

Then the multi-orthognoal polynomials of type II Plfj(x), [3J, [I], [IB], |1D] are polyno- 
mials of degree N — 1 such that 

P^ l {x)x j w {x)dx = 0, < j < N - 3, 

(4.9) 

PN,i( x ) w m(3c)dx = -27ri5i m , I, m = 1, 2. 

+ 

Remark 4.1. More accurately, these are in fact the multi- orthogonal polynomials with 
indices (N — 2, e*/) ; where t\ = (0, 1) and e 2 = (1, 0). 

We will now define the multi-orthogonal polynomials of type I. Let P^^x) be a function 
of the following form 

2 

P N,ii x ) = B N j(x)w + 2J (hix + A NA i) w k , (4.10) 

k=l 

where B^^x) is a polynomial of degree N — 3 and A^ t k,i is independent on x. Moreover, 
let Plj^x) satisfy 

/ P^(x)x j dx = 0, j = 0,...,N-l. 

Jr + 
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Then the polynomials Bnj(x) and 5kix + are multi-orthogonal polynomials of type 

I with indices (N, 2, 1) for I = 1 and (N, 1, 2) for I = 2. We will now show that Pffi and 
Pp/ 2 exist and are unique if both (Lat-i, Ln- 2 ) 1 and (L/v_3, Ljy_4)i are non-zero. 

Lemma 4.2. Suppose both (L/v-i, Ljq- 2 )i (Ln-3, Ln-a)i are non-zero, then the poly- 
nomials Pjfi and Ptf2 ex ^ s t an d are unique. 

Proof. Let us write Ljv_; as 

= 3- (qi(x)x(t - x)w) wT 1 + Ri(x), 1 = 1,..., A. 
ax 



for some polynomials of degree N — I — 2 and -R; of degree 1, then by Lemma 14.11 
we see that both Qn-i and Qn-3 are invertible and hence the composition Qn-iQns is 
also invertible. In particular, there exists a 2 x 2 invertible matrix with entries and 
polynomials q~i of degree N — 3 such that 

Ln-I = (qi(x)x(t - x)w) W' 1 + Ci- 2 ,\Ln-1 + Q-2,2^JV-2, / = 3, 4. 

Then by the first condition in (14.91) . we see that 

/ PNjWj(x)dx = {Pn\i, Q,iL/v-i + ci )2 L N ^ 2 ) 1 = -2m5i 



1 = 1,2, j = 1,2. (4.11) 



As the matrix with entries is invertible, we see that the linear equations (I4.1ip has a 
unique solution in the linear span of L/v-i and L/v-2 if and only if (Ljv-i, Ln_ 2 ) 1 ^ 0. □ 

As we shall see, existence and uniqueness of PjJ l would imply that the multi-orthogonal 
polynomials of type I also exist and are unique. As in the multi-orthogonal polynomial 
together with the skew orthogonal polynomials form the solution of a Riemann-Hilbert 
problem. Let Y(x) be the matrix 



Y(x) 



( 7Tjv,l(x) 
K7T Ar _2,l(x) 



C (n Ntl w ) C (n NA wi) C (n NA w 2 ) \ 



V 



Dll 

TV, 2 



X 



(4.12) 



where k is the constant 



2m 



Am' 



K 



Mh 



JV-1,1 



h 2 j-l,l — (7T2_j-2,l) 7T2jf— l,l) i • 



(4.13) 



and C(/) is the Cauchy transform 



lis) 



ds. 



x 



(4.14) 
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Then by using the orthogonality conditions of skew orthogonal polynomials and multi- 
orthogonal polynomials, together with the jump discontinuity of the Cauchy transform, 
one can check that Y(x) satisfies the following Riemann-Hilbert problem. 



1. Y(z) is analytic in C \ M. + , 

(\ w wi(z) w 2 (z)\ 

10 

1 

\0 1 / 



2. Y + (z) = Y_(z) 



3. Y(z) = (l + 0(z- 1 )) 



fz N 



\ 



y-N + 2 



z e 



\ 



7 



(4.15) 



z — > oo. 



Similarly, multi-orthogonal polynomials of type I can also be arranged to satisfy a Riemann- 
Hilbert problem. Let e be the operator 



[x-y)f(y)dy. 



(4.16) 



First note that the functions ipj(x) = e(7Tjitu) for j > 2 can be express in the form of 
(14.101) . By Lemma [4.11 we can write 71^1 (x) as 



Kj,i( x ) — ^ (Bj(x)x(t - x)w) w 1 + A jA L N _ 3 + Aj, 2 L N -4, 

Then we have 

e(-K jt iw)w = Bj(x)w + Aj^Wi + A jt2 w 2 . 

where Bj(x) is a polynomial of degree j — 2. 

Let X{z) be the matrix value function defined by 

(-^fC^ N . 2W ) !fB N . 3 ^fA N ^ 2 \ 

-fc(M 



(4.17) 



X(z) 



M R 

—Bn 



M a 



M A 



(4.18) 



Then by using the orthogonality and the the jump discontinuity of the Cauchy transform, 
it is easy to check that X~ T (z) and Y(z) satisfy the same Riemann-Hilbert problem and 
hence the multi-orthogonal polynomials of type I also exist and are unique. 

We will now show that the kernel Si(x,y) given by (I2.8p can be expressed in terms of 
the matrix Y(z). 

Proposition 4.2. Suppose (L A r„ 3 , L A r_ 4 ) 1 (L/v-i, L N _ 2 ) 1 7^ and let the kernel S\(x,y) be 



N-l 



Si(x,y) - - ^ r j (x)w(x)/j ljk e(r k w)(y), 

j,k=0 



(4.19) 
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where rj(x) is an arbitrary degree j monic polynomial for j < N — 2 and rj(x) = ivj^x) 
for j > N — 2. The matrix fi with entries fij k is the inverse of the matrix D whose entries 
are given by 

(P) jk = (r j ,r k ) 1 , j,k = 0,...,N-l. (4.20) 
Then the kernel Si(x,y) exists and is equal to 



/A 

o 



w 



(4.21) 



Proof. First note that, since (L^-i, L N _ 2 )i ^ and (L N _ 3 , L N _ i ) 1 ^ 0, the skew orthog- 
onal polynomials 7Tjv-z,i exist for / = —1, . . . , 2. In particular, the moment matrix D with 
entries 

(fi) k = {xi,y k ) x , j,k = 0,...,N-l 

is invertible. Since the polynomials ^N-j,i exist for j = —1, . . . , 2, the sequence r k (x) and 
x k are related by an invertible transformation. Therefore the matrix D in (14.20p is also 
invertible. As the matrix D is of the form 

©o \ / 1 

hN-^JJ' J \-l 

where Do has entries (rj,r k ) 1 for j, k from to N — 3. From this, we see that the matrix 
fijk is of the form 

As in [28], let us now expand the functions xrj{x) and xefaw). 



N-l 



k=0 

N-l p / , „> P i 

xe 



d jk e(r k w) + d jtN ip N + d j>N+1 — ^— + d j>N+2 - 



(4.23) 



A _ (j 



Then the coefficients Cjfc and djfc for j, k = 0, . . . , N — 1 are given by 

N-l N-l 



1=0 1=0 



Therefore if we let C be the matrix with entries Cj k , j, k = 0, . . . N — 1 and X> be the matrix 
with entries dj k for j, k = 0, . . . N — 1, then we have 



C = D x /i, £> = (/iD!) T , (4.24) 
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where Di is the matrix with entries (xrj, r k ) 1 for j, k — 0, . . . , N — 1. From (I4.23p . we 
obtain 

N-l 

(y - x)S 1 (x,y) = ^■K N)l (x)w{x)^ N ^ ld e(r j w) (y) 

j=0 

- y rj{x)w{x)^ jk dk,NVN{y) + d k ,N+i — r~r + d k.N+2 — r~r 
jj^o V w ^ y) ^ w ^ y) ) 

+ r T (x)w(x) (C T ji — \xD\ e(rw)(y), 
where r(x) is the column vector with components r k (x). Now by (I4.24p . we see that 



j 



which are equal as // = — fi. Now from the form of /i in (I4.22p . we see that ///v-i. 
^•,jv-2^jv_i i- Hence we have 

N-l 

^ 7!- NA (x)w(x)fi N - ltj e (rjw) (y) = h N 1 _ ll Tx N:1 (x)w(x)i) N - 2 (y). (4.26) 

3=0 

Let us now consider the second term in (I4.25p . As in (14.171) we can write e(r k w)w as 

e{r k w)w = q k (x)w + D ktl Wi + D kt2 w 2 , 

where q k (x) is of degree k — 2. Then from the form of in (14.101) and the orthogonality 
condition (I4.9p . we see that the coefficients d ktN+h I — 1, 2 are given by 



4,JV+i = Am = -— / P N I i(x)e(r k w)wdx. 



For A; 7^ N — 1, the polynomial g/c is of degree less than or equal to N — A, while for 
k = N — 1, Bn in (14.171) is a polynomial of degree N — 2 and hence the coefficient d kt N is 
zero unless k = N — 1. For = iV — 1, it is given by the leading coefficient of qN-i divided 
by the leading coefficient of B N . Since 

7rjv-i,i(x) = ^ (gAr_i(x)x(t - a;)™) w~ l + D k:1 L N ^ 3 + D fej2 Ljv-4, 

we see that both the leading coefficient of g7v_i(x) and B N is -g. Hence <i fe A r is S^n-i- 
This gives us 

jV-l 

} u rj{x)w(x)^ jk d kyN ip N (y) = -h N 1 _ 11 TX N ^i(x)w(x)%l) N (y). 

k,j=0 
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To express the second term in (I4.25P in terms of the multi-orthogonal polynomials, let us 
now express P^ in terms of the polynomials r k . Let us write P^ = J2f=o a j r j( x )- Then 
we have 



N-l N-l „ 

P"i{x)e{r k w)wdx = V" a 3 - (D x ) ifc , V] / PN,i(x)e(r k w)wdx^i kj a r 



j = fe=0 



Hence P^^x) can be written as 



jv— 1 / /. \ jv— 1 

^0*0 = ^ f / PN,l(x)e(r k w)wdx/I k j) r i( a; ) = d k,N+l^kfj{x) 
k,j=0 ' k,j=0 



Therefore the second term in (I4.25P is given by 

N - 1 / Pni(v) P N2(y) 



j,k=0 



w(y) ' w(y) 

2 

2tt/ 



! 2 

■hN-i,i*N-2A x )' , l>if(y) + ^z2 p N,i( x ) w ( x ) wl (y) p N,i(y)- 



i=i 



From this and 04.261) . we obtain 



(j/ - x)Si(x,y) = /l Ar 1 _ 1>1 7T A r i i(x)w(x)V'Ar-2(y) + /i JV ] _ lil 7rjV-2,l(^)w(^)V'Jv(2/) 



By using the fact that Y l (y) = X T (y) and the expressions of the matrix Y (I4.12p and X 
(14.181) . together with (14.131) . we see that this is the same as ( 14.2 ip . □ 



4.3 The kernel in terms of Laguerre polynomials 

We will now use a result in j9] to further simplify the expression of the kernel S\(x,y) 
so that its asymptotics can be computed using the asymptotics of Laguerre polynomials. 
Let us recall the set up in [5]. First let Y(x) be a matrix satisfying the Riemann-Hilbert 
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problem 



1. Y(z) is analytic in i, \ ir\ + , 

/l wo{z) w\{z) ■■■ w r (z)\ 



2. Y + (z) = YJz) 







o 
i 



z e 



3. Y{z) = (I + 0{z- 1 )) 



-n+r 



V 



z — > oo. 



V 



and let /Ci(x, y) be the kernel given by 



27ri(x — y) 



(0 w (y) ... w r (y))Y- 1 (y)Y + (x) 







(4.27) 



Let ^-,2(3;) be the monic orthogonal polynomials with respect to the weight wo(x). Let 
K,q(x, y) be the following kernel. 



JC (x,y) = w (x 



(4.28) 



h n -i,2( x - y) 

where hj ;2 = J °° 7T^ 2 Wodx. Let 7r(z), v(z) and u(z) be the following vectors 

7r(z) = (7T n _ r)2 , . . . ,7T„_i )2 ) , v(z) = [Wi, . . . , W r ) W Q , u(z) = [I - /C Q ) V (z) (4.29) 



and let W be the matrix W = J R 7r(z)f T (z)u7o(-2 ; )^- Then the kernel JC\(x,y) can be 
express as [9] 

Proposition 4.3. Suppose j R+ p(x)wi(x)dx converges for any polynomial p(x) . Then the 
kernel JCi(x,y) defined by ffJFfy is given by 



fCi(x,y) - K, (x,y) = w (x)u T (y)W l ir(x). 



(4.30) 



Remark 4.2. Although in JM/, the theorem is stated with the jump of Y(x) on R instead 
ofK + , while the weights are of the special form w = e ~ NV ( x \ wj(x) = e ajX , where V(x) is 
an even degree polynomial, the proof in JM/ in fact remain valid as long as integrals of the 
form j R+ p(x)wi(x)dx converges for any polynomial p(x) . This is true in our case. 
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We can now apply Proposition 14.31 to our case. In our case, the vectors tt(z) and v(z) 
are given by 

7r(z) = (LjV-2 ^7v-i) T , v(,z) = (wi w 2 ) T Wq 1 , 



while the matrix W is given by 

W = 



(Ljv_2, Ljv-3/i (Ln-2, -^JV-4/i 
(Ln-1, Ljv-3/i (Ln-1, Lpf-4/i 

By corollary 14.11 we see that (L N _ 2 , Lns)i = and hence the determinant of W is 

det >V = (LjV-2, £/V-4)i (£/V-l, ^JV-3)i • 

Suppose det W 7^ and that the multi-orthogonal polynomials exist. Let us now 
consider the vector By the Christoffel-Darboux formula, we have 

N-l j , \ 

u(z) = (/ - k-0 „(*) = - ir^ (^i^^J , 

where /i^o = (Lj,Lj) 2 . We shall show that L7V-3 and Lat-4 can be written in the following 
form 

L N _ t = —(qi(x)x(t- x)w)w~ l + Qi-2^n+i,i + Qi-2,2ir N ,i, Z = 3,4, (4.31) 

for some polynomial qi(x) of degree N—l. By Lemma |4TT| we see that if (L N _ 3 , L N _ 4 ) 1 ^ 0, 
then the map Qn-3 in Definition 14. II is invertible. Therefore if the restriction of the map / 
in Definition 14.11 is also invertible on the span of 7Tjv+i,i and 7Tjv,i, we will be able to write 
L7V-3 and L^-a in the form of (I4.3ip . 

Lemma 4.3. Let f be the map in Definition \4-l\ and let g n be its restriction to the span 
of Vjv '+1,1 and ttn,i- Then g n is invertible. 

Proof. Suppose there exist at and a 2 such that 

o-ittn+i 1 + ct2^N 1 = ~r~ (q(x)x(t — x)w) w^ 1 
ax 

for some polynomial q(x) of degree at most N—l. Then we have 

(x J , aiTTAr+1,1 + a 2 vr A r i i) 1 = {x 3 q(x)) 2 = 0, j = 0, . . . , N - 1. 
As the degree of q(x) is at most N — l, this is only possible if q(x) = 0. □ 
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The composition of Qns and g' 1 will therefore give us a representation of L^-z and 
Lat_ 4 in the form of (14.311) . By using this representation and the fact that qi(x) is of degree 
at most N — 1, we see that 

N-i r , x 

^ fin n 

= V] -7 {{Lj, qi) 2 + (Lj, Ql-2,l^N+l,l + (5/-2,27TjV,l)l) ) 

fin n 

i=o J ' U 
= qi(x). 

Therefore the vector u(x) is given by 

u(x) = w(x)wq 1 (x)Qe (n N+1< iw n Ni iw) T 
where Q is the matrix with entries Qij- We will now determine the constants Qij- 

Lemma 4.4. Let 1x^,1 and 7Tjv+i,i be the monic skew orthogonal polynomial with respect to 
the weight w(x) and choose 7Tjv+i,i so that the constant c in \4-8\) ^ s zero. Then the vector 
u(y) in given by 

u{x) = w(x)w l {x)Qe (7Tjv+i,iw kn,iw) T , (4.32) 

where Q is the matrix whose entries Qij are given by 

M (L N _i, L A r_j_ 2 ) 1 



Qi, 



i.i 



2/ijv-i,o 



Qv = (Mt -(N+ M)) (LAf - 1 ; Ljv — 2)1 - M {Ln - 2 ; Ln - 1 - 2)1 

2/liV-l,Q 2/ljv_20 



(4.33) 



Proof. First let us compute the leading order coefficients of the polynomial qi(x) in (14.311) . 

Let qi(x) = q^N-i^^ 1 + qi,N-2X N ~ 2 + 0(x N ~ 3 ), then we have 

(qi(x)x(t - x)w) w' 1 = y%n-iX N+1 , 
From (I4.3ip . we see that 

%N-l = -JfQl-2,1- ( 4 - 34 ) 

On the other hand, by orthogonality, we have 

(Ln-i, Ln-i)i = (Ljv-i,9i) 2 — %,jv-i^jv-i,o- 
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Therefore Qi-2,1 is given by 

M (Ljv-i, Ln-i) 1 

Vt-2,1 - ^7 

Let us now compute (5i,2- By taking the skew product, we have 

(Ln-2, Ln-^ = qi : N-i (Ln-2, % x ) 2 + %,Ar-2^Jv-2,o 
Now from (I4.4p . we obtain 



(M-N + j + 1)(j + 1) 
M 



(4.35) 



Hence qi,N-2 is equal to 



(L jy _ 2 ,L jy - i ) 1 (M-lXjV-1) 

9/.JV-2 r Qi,n-i 77 (4.db) 

riN-2,0 m 



By using the expansion (14.81) of the skew orthogonal polynomials in terms of L^-k, we 
have 

A, 17 , f , Mt M \ , 

+ I ~ m )Qi,n-i Y qi ' N ~ l + T qi ' N ~ 2 ) '° 

By substituting (14.361) into this, we obtain 

Ql _ %2 = (A (M + N) + —) - 

This proves the lemma. □ 

Note that the matrix det(<5W _1 ) = M 2 /(4/i7v_ 10 /i7v-2,o) and hence QW~ l is invertible. 
From Lemma [4.41 and (I4.30p . we obtain Theorem 12. 2[ 

We will now further simplify the expression for the correlation kernel S\(x, y). First we 
make the following observation. From the definition (I4.19P of the kernel Si, it is easy to 
check that the following identity is true. 

-|U(*,y) = ^Si(y,*). (4.37) 
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We can use this simple observation in the expression fl2.10p that we obtained in the last 
section. First let us write the Christoffel Darboux kernel terms of the solution of a Riemann- 
Hilbert problem, we have 

= (SS) ' mhv> (0 *> z " fe)z + w • (4 - 38) 

where Z(x) is the following matrix 

z <*> = (« L f (,) k c cfT o) » i) < 4 - 39 > 

\K0,N-1-LN-i[X) ^O.TV-lV \J-'N-1Wq) J 



where <7 3 is the Pauli matrix cr 3 = diag(l, — 1), K 0n = — and C(/) is the Cauchy 
transform. 

It is well-known that if the logarithmic derivative of Wq is rational, then the matrix 
Z(x) in (I4.39P satisfies a linear system of ODE 

d 

— Z(x) = A(x)Z(x) (4.40) 
ox 

for some rational function A(x). By using (I4.4p and the recurrence relation of the L n 
n ( M-N-\-Mx\ n(n-l) ( M — N — 1\ 

L - = -m 2 + ; - - W 1 + — s~ (4 - 41) 



we obtain the matrix A as 



_M , M+N N_ 



x 2 2x 



Let us now consider the derivative 9K2 Qy' y ^ . We have, from ( I4.38P 
dK 2 (x,y) t-2y 



d V (xy(t - x){t - y)Y 4m(x - y) \0 



Z- + \y)Z + {x) fl 
I 1 ' Am(x - y) V° 



; v(t-y) y , Z^MZ+Qe) A\ Z+^y^Z+Cx 
1 ~^ ^' 11 ' 2m(x-y) 2 [OJ { ' 2m(x - y) 



(4.43) 



Similarly, ^K 2 (y,x) is given by 
dK 2 (y,x) t-2x _ 



dx 



Z-\x)Z + {y) fl\ 
' Am(y-x) V / 



(xy(t - x)(t - y)Y 
f x(t-x) y ( Z-\x)Z + {y) (l\ Z-\x)A(x)Z + (y) f\ 

+ \y(t-y)J l lU > 2m(y-xy \0j + lU > 2m{x - y) V° 



(4.44) 
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The sum of these two derivatives is therefore given by 

h~ 1 N _ 1 w{x)w(y)\ - L N (x)L N (y)M 



dK 2 (y,x) dK 2 (x,y) , 



dx dy 

IT I \ T I \ T I \T l^f M l \ Mt+(M + N) 

+ (L N ^(y)L N (x) + L N ^(x)L N (y)) —{x + y) - 



2 v 2 



- LAr_i(x)LAr_i(|/)Ar 



By using the recurrence relations of the Laguerre polynomials, we obtain 
8K 2 (y,x) , dK 2 {x,y) x J (N-1)(M-1) T 

— di — — dy — = h o,N-i w ( x Mv) I 2M l k-zA x > y) 

M T , . (M + N) — Mt . , 

+ Y L N+i,N-i{x,y) H L N _ hN (x,y) 

where L n , m (x,y) is given by 

L n , m {x,y) = L n (x)L m (y) + L m (x)L n (y). 
Let Ki(x,y) be the correction kernel in fl2.10p . 

Ki(x, y)=e {ii N +i,iw tt n ,iw) (y) M mAn±m) ir(x)w(x) 



Then we have 



dKt(y,x) , dK x (x,y) ,,,J M 

j< + j< = w(x)w(y) - — L NjN ^ 2 {x,y) 

dx dy y 2h , N - 2 

M Mt — (N + M) 
- 777 ^JV+i,Ar-i(^,2/J H 777 L^N-^x^y) 



M (L N ,L N - 2 ) 1 M (L N+1 ,L N _ 1 ) 1 



2^0,^-2 (Ln-I, Ln-2/i 2/i 0i jv_l (I/jv-1, Ljv_2/i 



( M (L N+1 ,L N ^ 2 ) 1 Mt-{N + M) (L N ,L N - 2 ) 1 , 
+ I 77 7 v 7 77 7 r ) L N -i,N-i(x,y) 

Then by using ' 9l5 ' 1 ^' :r - > + 9Sl Qy' v ^ = and the formula for ft ,n 

_n!(n + M-A^)! 
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we obtain the following. 



it r v _ (iV-l)(M-l) 

\lvjV+l, J^N-lli — TT2 WAT, > 

/r r v _ Mt-(iV + M) 
In particular, the correction term K\ (x, y) in the kernel can be written as 



(4.50) 



M 



2h 



0,N-2 



e (L N w) (y)L N _ 2 w(x) 



(L N , L N _ 2 ) l 
M 



(e (Ln^w) (y)L N _ 2 w(x) - e (L N „ 2 w) (y)L N _ 1 w(x)) 



2k 



0,N-1 



% / \ T f \ Mt — (N + M) 
e (L N+1 w) {y)L N _ 1 w{x) H — e {L N w) {y)L N _iw(x) 



2h 



0,N-1 



(4.51) 



5 Derivative of the partition function 

In this section we will derive a formula for the derivative of determinant of the matrix D 
given in (I4.20p . We have the following. 

Proposition 5.1. Let D be the matrix given by ^W^, where the sequence of monic poly- 
nomials rj(x) in lji4-20\ ) is chosen such that rj(x) are arbitrary degree j monic polynomials 
that are independent on t and rj(x) = nj,i(x) for j = N — 2, N — 1. Then the logarithmic 
derivative of det D with respect to t is given by 



9, , . 
-log det 



Si(x,x) 
t — x 



dx, 



(5.1) 



where Si(x,y) is the kernel given in l[4-19\ ). 

Proof. First let us differentiate the determinant det D with respect to t. We have 



d_ 

ot 



( &D, 



det D = det 



+ det 



D 



01 



D 2n 

-1,1 

/ D 00 d t D m 



\D 2n -ifi d t D 2n ^ii ■ ■ 
+ ■ ■ • + det j : 



A),7V-1 \ 

Dn-i,n-i ) 

Dqm-1 \ 



Dn-i,n-i ) 



N-l 
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Computing the individual determinants using the Laplace formula, we obtain 

d JV_1 

— det D = det D ^ ^A^i- ( 5 - 2 ) 

i,j=0 

As Tj(x) are independent on t for j < N — 2, the derivative dtDij is given by 

^--Kte''), + (-^)>-i(fe^-fe-> 1 

for i, j < N — 2. For either i or j equal to iV — 2 or N — 1, we have 

) n x I 1 f / n N-2,l \ I 7TjV-l,l 

d t A,JV-i = 0jv-2,i I - 2 ^\ t-g ' 7rjv ~ 1 ' 1 / ~ \ t - x ' 7FiV ~ 2 ' 1 

+ (d t 7T N - 2 ,l, 7TJV-I,l)i + (tTJV-2,1, ^^-1,1)1 J • 

Note that by orthogonality, the last two terms in the above expression are zero as dtitN-x,i 
is of degree N — 2 and d t 7i N _ 2jl is of degree N — 3. Applying the same argument to d t D i;N _ 2 , 
we obtain 

in X (I n N-2,l \ I KN-1,1 

Ot^>i,N-l = —0N-2,i - _ - ,71" AT- 1,1 1 - - _ - , TTjV-2,1 

Q n ^/V-l.i / / TTjV-1,1 \ / 7l N _2,l 
Ot^>i,N-2 = — ( — ,^N-2,1 ) - { ,^N-1.1 



As D is anti-symmetric, the derivatives dtD^-i^ and d t DN-2,i are given by dtD^^i. 
—d t D i N _i and d t D N _ 2ji = —d t D iN _ 2 . From these and (15. 2p . we obtain 



a / AT-3 

a 



-detD = detD 

• lJ= ° X 71 x (5.3) 

,, I / \ / 7rjv_i i \ \ 

-! - ( \ t _ % i n N-l,l j - ^ ~ _ ~ ) TTjV-2,1 ; J /'.V-2..V-J 

where we have used the anti-symmetry of D and \x to obtain the last term. From the 
structure of the matrix \i in (I4.22p . we see that the last term in (15. 3p can be written as 

2 ( ( JT^'^- 1 '^ ~ i^fzf^ N - 2 ^ ) i ) VN-2,N-1 
N-l JV-1 , , N-l N-l I i 

E E (r^' r v ^+ E E\^' r v a*« 

i=0 j=N-2 X ' 1 i=JV-2 j=0 X ' 1 

From this, ( 15. 3 p and the expression of the kernel in (14.191) . we obtain ( 15. ip . □ 
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6 Asymptotic analysis 



We will now use the asymptotics of the Laguerre polynomials to obtain asymptotic ex- 
pression for the kernel Si(x,y). We will demonstrate the asymptotic analysis with the 
assumption that M/N —7-7 = 1. The same analysis also be applied for general 7. 

Let us recall the asymptotics of the Laguerre polynomials in different regions of R+. 
These asymptotics can be found in [70]. When 7 = 1, the end points b- and b + of the 
support of p in (11.11) are and 4 respectively. 

Let n = N — j , where j = 0(1). Let us define the function tp(x) to be 

where the contour of integration is chosen such that it does not intersect the interval 
(—00,4). We can now use <p(x) to define two maps /„ and /„ in the neighborhoods of 
and 4. 

Definition 6.1. Let (p+(x) be the branch of <p(x) in the upper half plane and let 5 > be 

sufficiently small, then the maps f n is defined in a disc of radius 5 center at 6_ in as follows. 

f n =^ + {x)-v + m)\ 7 = 1- ( 6 -2) 

In the above expressions, the map f n is defined to be the analytic continuation of the 
expression in the right hand side to the whole disc. 

Similarly, the map f n in a neighborhood of 4 is defined to be 

2 

fn = (|"¥> + (*)) 3 ( 6 ' 3 ) 

As in JTm , these maps are conformal inside the small discs around b± >n , provided 5 is 
sufficiently small. In particular, they behave as follows as near b± >n . (Recall that b- t7l = 
and b +tn = 4 when 7 = 1.) 

fn = -n 2 x~f n (x), f n = (J) 1 (x - 4)/ n (x), (6-4) 

where the maps f n (x) and f n {x) are of the form 

} n {x) = 1 + 0(x), f n (x) = l + 0(x- 4). 

Note that f n (x) and f n (x) are independent on n. 

We can now state the asymptotics of the Laguerre polynomials. These can be found, for 
example, in [70]. When M — N = a = 0(1), the Laguerre polynomials have the following 
asymptotic expressions in different regions on R + (See [70]). 
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Proposition 6.1. Suppose M — N = a = 0(1), then for n = N — j and j = 0(1), 
there exists sufficiently small 5 > and the Laguerre polynomials L n (x) in i4-4\) have the 
following asymptotic behavior on R + as n, N, M — > oo. 

1. Uniformly for x G (0, S], 



n w /n Us 2y^F(-l) n (-/ n )3 (m n _!) » . 

L « ( x I ( ( — T x ) ) = — -r I sin (J a ( 2(-/ n 

%i (4 — x) * 



N J \ \N 



~ , i 

2 



n w /n \\§ 2 (m n _i) g / / 7T 



x 



(1 + (AT- 1 )) + cos + inp+ - ^) O^ 1 ) j 



(6.5) 



x (1 + 0(N- 1 )) + cos (2(-/„)i) (1 + 0(N- 1 )) 

where J a is the Bessel function and £ is the following 

1 (x \ ira 

C = -(a + l)arccos(^--lj -— . (6.6) 



T7ie branch cut of arccos(x) m £/ie above is chosen to be (— oo, — 1] U [l,oo). The 
constant K n is given by K n = 

2. Uniformly for x G [5, 1 — 5], 



(6.7) 



where rj± are given by 

V± = "*■) arccos ^ ~~ lj (6-8) 

5. Uniformly for x G [4 — 5, 4 + 5] , 

x (1 + OiN- 1 )) - \f n \~* S in V+ (x)Ai'(f n )(l + OiN- 1 )) 



(6.9) 



Remark 6.1. Reader may notice that the asymptotic formula presented in Proposition \6.1\ 

is different from the ones in Yffl$. This is because the weight in |77¥ that is relevant to us 



is Wq(x) = x a e 4nx and a rescaling of the variable from x to y = jj-x is needed to obtain 
the formula in Proposition \6.1\ from the results in 



We will now use the asymptotic formulae in Proposition 16.11 to compute the skew 
products of the form (L N _j, L N _ k ) v We shall follow the ideas in |3T], [33] and 
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6.1 Asymptotics for the skew products 

Throughout the analysis, we will assume that t is of finite distance from [0,4). We will 

2 

also need to consider the case where \t — 4| is of order iV~3. We shall organize the results 
according to the different regions on R + . We will divide R + into four regions (0, N~*], 
[N~5, 4 - N~ £ ], [4 - N~ £ , 4 + N~ e ] and [4 + N~ £ , oo) for some e > 0. These regions are 
called the Bessel region, the bulk region, the Airy region and the exponential region. 

6.1.1 The Bessel Region 

Let us define L n (x) to be 



L n (x) = L n \ W [j7 x J (««n-l) 2 



K ' a7 x ) (S GH) 2 (* - ^ 2 



(6.10) 



From (16. 6p . we see that as x — > 0, sin£ and cos£ have the following behavior. 

a 4- 1 

sinC = l + 0(x), cosC = — ^— v^(l (6.11) 

Let us now express the integrals in the Bessel region in terms of the variable f n . From 
(16. 2p . we have 

™ _ _ln_ |i L^-lD I ^ 

nA Yin 2 + n 4 



(6.12) 

Vat V AT 77 I tn 2 ) I In 2 

As x < iV~5 in the Bessel region, we see that the term /^/n 4 is of order 0{n~ l ) at most. 
Then from the asymptotic formula (16.51) . we have the following asymptotic formula for 
L n (x) in the Bessel region. 



- (-l)-V^ni ( ■>« (2(-/-) i ) / t 1|0 ~ 

" + 1 V (of t \l\ ( 1 , fc 2>° i r fc 2j/2 



n 



2j 



(6.13) 



n 2j 

3=1 



for some constants fcij and that are bounded in n. By using the asymptotic formula 
for the Bessel function in PQ, we obtain the following. 
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Lemma 6.1. Let k>Os>v>0 and a > 0, then as s — >■ oo, we have 

J J a (u)u k du = 0(1) + O (s fc -^ , J J> a ( u )u k du = 0(1) + (s fc ~^) . (6.14) 

Let si = (s) as s — » oo ; and let B\ = J a , B 2 = J' a , then as s — > oo, we have the following 
for the double integrals. 

[ Bi(u)u k f 1 B j (v)v l dvdu = 0(s k+l ) , i,j = 1,2, k + I > 0. (6.15) 

JO Ju 

Proof. Integrating the first formula by parts, we obtain 

J a (u)u k du = ( - [ J a (v)dvu k ) +k [ ([ J a (v)dv J u k ^du 



V \J u 



By (9.2.1) and (11.4.17) of [T], we see that ff J a (v)dv = 0(s 2) as u — > s and s — > 00. 
The function J s J a (v)dv is therefore bounded and behaves as 0(u~^) as u — > s. Therefore 
we have 

jf J a (u)u k du = k (j J a (v)dv^ju k ~ 1 du + 0(l) + 0(s k ^) = 0(l) + 0(s k -^). 

This proves the first equation in (16. 14p . Integrating the second equation by parts and using 
Ja(0) = for a > 0, we obtain 

J' a (u)u k du = J a (s)s k - J a {v)v k - k j J a (u)u k ~ l du. 

J v 

By using the following estimates for any C > 0, 

sup \y~ a J a (y)\ < 00, sup \y/yJ a (y)\ < 00, 
velo, C) w6[ c,oo) 

sup |y a+ J a (y)| < 00, sup |Vj/J a (j/)| < 00, 
j/e[o,C) ye[C,oo) 

we obtain the second equation in (I6.14p . The estimate (I6.15P now follows immediately from 
(EUD and f l06]) . □ 

We can now compute the integrals in the skew products in the Bessel region. 
Proposition 6.2. The single integral involving L n (x)w(x) in the Bessel region is given by 



_ 1 

/* N ^ / -1 \ n /o 

1* / L n (x)w(x)dx = l ln = l ~ ' \1 7r +0(n~l). (6-17) 
The double integral in the Bessel region is given by 

iy/K n -iK m -i / L n (x)w(x) / L m (y)w(y)dxdy 

Jo Jx 



U/Cn-1 



(-l) m +"27T ^. . 

nr 



(6.18) 
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Proof. Let us first prove the single integral. We have, by (I6.10p . 



i 



-A — 1 v^V 



(m n _i) 2 / L n (x)w(x)dx = I L n [—xjdx= — \ L n (x)dx. 



N \ . n 



From (I6.13p . we have 



f^^(r(«(-^l , - i,v ' 



(2n) 2 J 



a + 1 t/ / \ / u k 2 au ^ (-iyk 2 ,jU 23 +1 \ \ 

I — J n {u) — t^+> „. ' dw 

4 aV ; n 2 ^ 2 2 -?n 2 ^ +1 y J 

where we have changed the integration variable from x to u = 2{—f n )^ and u + is the 



corresponding upper limit u + = 2^j—f n (^f^ = O ( n *)- By Lemma |6TT| we have 

jf " L n (x)dx = ^-^=^ f jf + J a (w)dw + O (V§) 

By using the asymptotic formula for the Bessel function (pQ, (9.2.5), (9.2.9) and (9.2.10)), 
we have obtain the following estimate 

I J a {u)du = O (u^) (6.19) 

J u+ 



By using this and the fact that J °° J a {u)du = 1, we arrive at (I6.17p . 

Let us now compute the double integral. Let n = N — j and m = N — k, where j 
and k are finite. As in the single integrals, we will change the integration variables into 

u = 2\J —f n and v = 2\J—f m . Let v+ = 2\J —f m (\/N /n) = O (mi^ be the upper limit in 
the variable v and v{u) be the value of v when y = —x. Then by using (I6.12p . we have 

^yH i+ ^ +o (3) (6 - 20) 

Now the double integral is given by 

iy/Kn^Km^ / L n (x)w(x) I L m (y)w(y)dxdy = — / L n (x) / L m (y)dydx. 

x m x 

By using Lemma 16.11 and ( I6.16P . We see that it is of the order 

r^r . r 3 ^ . (-l) m+n 2n f u+ f v+ ( 

/ L n (x) L m (y)dydx = y ' / J a (u) J a (v)dvdu + 0(N-± 

Jo J^x V™nt JO Ju(u) V 
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To compute the first term, note that, since 

v'{u)J a {v{u)) ■■ 



d 

du 



J a (v)dv 



u(u) 



we have, by (I6.20p and mean value theorem, the following 

(in — n)u 



(l + Oin- 1 )) [Uu) + J' a (Q 



2n 



du \Ju(u) 



where £ M is between u and (u) is bounded, we have 



J a {u) / J a (v)dvdu 

(I Ju(u) 



2n 



\ 2 / r v + ^ 2 

J a {v)dv - / J a {v)dv 

V+ 



J'aitu) [ (m ^ +o(n-i))+0 (n- 1 ) ) / J a (v)dvdu. 



v(u) 



By using Lemma [6.11 and (I6.19p . we obtain 

/ J a {v)dv = O (n-t) , / J a (u)du = l + O (n~i ) , 

^ ' Jo ^ ' 

r i^r 1 + ° h)) + ° ^) r = ° 

This gives us 



_ j. 
n 4 



r u+ r+ i / _i 

/ ■/«(«) / Ja(v)dt;du = - + I n 4 

JO 7i/(te) ^ V 



This proves the proposition. 



□ 



This concludes the analysis in the Bessel region, we will now consider the integrals in 
the bulk region. 



6.1.2 The Bulk region 

First recall the following matching formula in the region 
Lemma 6.2. (Proposition 4-4 °f [33]) 



2n ' ' 



from 133 1 . 



1. Uniformly for x G 



VE $ 

2n ' U 



, as n — > oo, 



-f n j ( sinCJa ( 2 ( -/ n ) ) + cosCJ; ( 2 ( -/„ 

f-1 



(6.21) 



7T 



(cos F n (x) + L n sin F n (x)) + 0(n ) 



where i n 



16(-/„)2 



- and F n = r/ + + inip + — f . 
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2. Uniformly for x E [4 — e, ^c_] , as n — >■ oo, we /icroe 



cosr7+|/ n |4v4z(/„) - sin77 + |/ n | = -j=cosF n + ( — — - — - 

V 71 " Vn(4-a;) 



(6.22) 



This shows that throughout the bulk region, the function L n (x) in f)6.10p has the fol- 
lowing asymptotic behavior. 



L n (x) 



2VN 



- {cos F n (x)+E) 



(6.23) 



uniformly for x G 



VN N 
2n ' n 



y/nx*(A — x)a [t— ^x)' 

. , where E is an error term of order O fn _1 (4 — x)~^\ + 
O yn~ 1 x~sj , From this we can now compute the single integral. 

Proposition 6.3. Let c_ = 4 — N~ e such that e < ^. Suppose the distance from the point 



t to the interval 



VN N 
n ' n 



is greater than CN 6 for some C > and the distance from 



t to is finite. Then the single integral in the bulk region is of the following order. 



I i L n (x)w(x)dx = X 2 , n = O (n 



(6.24) 



Proof. By (I6.10p . we have 



K^i) 2 / i L n (x)w(x)dx 



n / « 



N Un 



L n (x)dx. 



From this and (I6.23p . we obtain 



cosF n (x) + E 



L n (x)dx = 2\ — I 
m \ N JvE X 3 Z {4 - x)l (t - §x) 



-dx. 



(6.25) 



Since F n = r} + + in<p + — j, we see that its derivative is given by 
FL ■■ 



(a + 1) in 

— , H 

2^x(A - x) 2 



x - 4 



a + 1 + n(A - x) 



2yJx(A - x) 



(6.26) 



where the plus subscript indicates the boundary value as y/x — 4/ yfx approaches the real 
axis in the upper half plane. 

This gives us the following estimate of the order of 1/F' n . 



F' 



2y/i 

ny/A — x 



1 + 



n(A — x) 



(6.27) 
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We can now integrate the first term in (I6.25P by parts to obtain 



cosF n (a:) 



— dx 



sin F n 



^ xi{A-x)\{t-%x)* F>f(4-x)i (t- §x) h 
d I I 



smF n (x) , 
M dx I ^1(4- x ) \ {t-^xY 



dx 



As the distance from t to the bulk region is at least of order O (N~ e ), the first term in the 
above equation is of order O ^n^ _1 j + O (n~i j, which is of order O ^n^ij as e < 1/20. 
By repeating the integration by parts procedure to the second term, one can verify that it 
is of order at most O (n 1 ^' 2 ^. This gn 



nves 



cosF n (x) 



v .V 



art (4 - 2)3 (t - ^x) 



—dx — O in 



Since the error term E is of order O (n 1 (4 — x) 2 \ +0 (n 1 x a j , we obtain the following 
estimate for contribution from E. 



n / « 



^-'^ xt(4-x)i(t-fx) 
This concludes the proof of the proposition. 



r dx = O (n * 1 



□ 



In particular, from the proof of Proposition 16. 3[ we see that 



cos F„ (x) 



y xi(A — x)i (t — jjx) 



—dx = O [n 



y e 



N N 



-C- 



n n 



(6.28) 



Let us now proceed to compute the double integrals. To begin with, we have the following 
lemma that will help us to simplify the results. (Proposition 5.12 in [33J) 



4-s 



Lemma 6.3. The function 9{x) defined by 

9{x) = - 

satisfies the following 



ds = —i(p+, x G [0, 4] 



9(x) — x9'(x) 



1 x 1 
arc-cos ( — — 1 



(6.29) 



(6.30) 
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Proof. By differentiating 6{x) twice we see that 

xO"(x) = ■ 

y x(4 — a;) 

As arccos(:c) is given by 

arccosla;) = / — = 

we see that 



arccos(x) = / ds, (6.31) 

A VI - s 2 



(9(x) — xO'ix)) = — — arccos ( l\ 

ax \2 J 



dx 

Integrating the above, we obtain 



(x) — x9'(x) = — arccos ^— — lj + C. 



For some integration constant C. By evaluating the above equation at x — 4, we see that 
C — 0. This proves the lemma. □ 

The lemma implies the following. (Proposition 5.13, [33J) 



Lemma 6.4. Let n = N — j , m = N — k where j and k are finite integers, then for 
x G 



n ' n 



, we have 

F m (x—^j - F n (x) = (m-n) arccos (| - l) + O (mi" 1 ) . (6.32) 

Proof. As in |33J, let us write the left hand side of ( I6.32p as 

F m [x-^j - F n {x) = [F m (x-^j - F m {x)) + (F m (x) - F n (x)) 

By repeat application of the mean value theorem, we see that the first term on the right 
hand side is given by 

„ / n\ . . x(n — m) „. . , x 2 (n — m) 2 .. ,. , 
F m (x-) - F m {x) = -i >-F' m {x) + \ > F" 



for some £ G [i,in/m]. From (16.261) . we see that for x G 



N_ c 

n ' n 



, we have 



y -F^(x) = (m - 7i)x0'(x) + O ' 



m "' \17ly/A — X 

x 2 (m — n) 2 // j 1 \ / 1 



(0 = o — — -r +o 



2m 2 x \m 2 (A-x)2j \m^/A - x 
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As 4 — x = N 6 and e < 1/20, we have 

F m far—) - F n (x) = (m - n)xfl'(x) + (F m (ar) - F n (x)) + O (m^ 1 ) 
= (m — n)x9'(x) + (n — m) 9(x) + O (m 5 " 1 ) 
= (m — n) arccos f — — l) + O (m2 _1 ) , 

where the last equality follows from lemma 16.31 □ 

We can now compute the double integrals in the bulk region. 
Proposition 6.4. The double integral inside the bulk region is of the following order. 



(6.33) 



(x)w(x) / L m (y)w(y)dxdy = J 2 

JN~2 Jx 

4 f ~ c ~ sin ((n — m) arccos (| — 1JJ 



f n c - sin ((n — m) arccos (§ — l)l f ^ 

/ U \u \ -dx + 0[n-i 

Jvn x(A -x)(t-x) V 



where c_ = 4 — iV e . 

Proof. The proof is similar to the computation in [33] and [3TJ. We have 



777, ?2 

^ = TP Ivn Ln ^ L L m{y)dxdy, 
Then by using (16. 23ft . we see that the double integral is given by the following. 

N N_ 

„ 4v / mn /"" c - cosF n (x) + /"™ c - cos F m (y) + E(y) 

J2 = — 77 — / r / — -dxdy 

N x i( A -x)\{t-%xY h* y*(4-y)4(t-ffy)* 

First let us compute the error terms. By changing the order of the integration, we have 

cosF n (x) fi c - E(y) 



Vn 3 
— X 



1(4-^)1 (t-%x)* J%* yl{A-y)l(t-^y) 
E(y) [nv cosF n (x) 



T dydx 

2 

T dxdy 



€ yl(4-y)i{t-% y yJ^ xl{A-x)\(t-%x) 

N 

= ^ r O fn-I) dy = (n-l) , 

% yl(A-y)\{t-% y y V 7 V 7 

where the last equality follows from (16. 28 p . The order of the other error term can be 
estimated similarly. Let us now consider the leading order term. Integrating by parts, we 
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obtain 



J: 



20 



cosF n (x) 



V^V 3 

— Xi 



cos F m (y) 



m\! /"« c - cosF ra (x) sinF m (gg) 



T dydx 



dx 



cos F„ (x) 



X 



f(4-x)3 (t-fx) 5 



sin F m (y) — 



+ O (n~i 



To evaluate the second term, note that for y in the interval 



n_ c 

m ' m 



we have 



smF m (y)- 



O 



5 5 
4 



d 2/ ^ 1/ i(4-y)4(t-ffiy)*/ V'"'- -!'^H 

Hence by interchanging the order of integration, we obtain 



cos F n (x) 



& xl 



(4-x)i {t-^xY h 



d 



sin F m (y) — 



O (m * j + O (m s 1+5 4 j = O (m * 



To compute the first term of J20, let us change the integration variable from x to —x. 
Then we obtain 



20 



m \ 4 / m 
n / /v^v 



cosF n sinF,, 



x 



^ F4(x)xf(4-x)i(4-ffix) i (t-2x) 



dx + 0(m 4 



Now by (I6.27P and the fact that both 4 — x and t — jjX are of order at least n £ , we obtain 
the follow estimate 



.*,=-2r-v 



cosF„ (^x) sinF m (x) 



dx + 0(m~ 2+3e ) + O ( m~i ) . 



\nJ J £EL mx(4 — x) (t — x) 

m 

To evaluate the integral, we use the angle addition formula for sine to obtain 



J, 



20 



n J \J^E mx(A — x) (t — x) 



my I j sin (F m (x) - F n (*x)) ^ | sin (F n (^x) + F m {x)) ^ 



y£E. mx(A — x) (t — x) 



+ 0(m 4 % 



48 



The first term can be simplified using (I6.32p while integration by parts shows that the 
second term is of order O (m~ij. This gives 

Cm c - sin (in — m) arccos (% — l)) 4 
J vw mx(A — x)[t — x) 

m 

As — c_ — — c_ is of order m -1 , we see that 

I™ ' ^ ((m - n) arccos (f-l)) ^ = Q +2 
,/iv c _ ma;(4 — x)[t — x) 

Similarly, we can change the lower limit to and add an error term of order O (m~ 2 ). 
This proves the proposition. □ 

As we are only going to consider the values of m and n with m — n = ±1 or ±2, we 
can simplify the double integrals further. First let us consider the case when m — n = ±1. 
A simple calculation shows that 

sin (arccos (| - l)) = vSEj? ( 6 . 34 ) 

From this we obtain 

sin (2 arccos (f - l)) = VgjEg^lj) . (6.35) 

We can now use these and residue calculation to compute the double integrals. 
Proposition 6.5. Let T\ be 

f 4 A sin (I arccos (f - l)) 

Ti= \ Ar \ A -t^ r^dx 6.36 

Jn c _ Nx(A-x)(t-x) v ; 

n 

Then the asymptotics of the double integrals in the bulk region are given by the followings. 

J 2 = ± ?1 T Fi + O (N-i) , n - m = ±1, 

NJt{t - 4) V / 

/ \ (6-37) 

Proof. We shall compute the integrals using Cauchy theorem. Since [y/x — 4) = i\/A — x 
on [0,4], we have, by using f)6.34p . the following 



in (arccos (f - 1)) f 4 1 

' ' dx = i / : dx, 



o x(A-x)(t-x) J 2 x /x{x-A) (t-x) 
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and similar relations for n — m = 2. The right hand side can be computed using Cauchy's 
theorem and we obtain 

n — m — 1; 



4 sin ((n — m) arccos (f — l)) 2y/t{t-4) 

dx - ; 



x{A-x)(t-x) f ( —7== — 1 I , n-m = 2; 



By fl634j) and (jOS) . we see that 



4 sin (/ arccos ff — l)) / 6 \ 



/ A^x(4 - ar) (i - x) 
From this and (I6.33p . (16. 34[) and (I6.35p . we have 



sin ((n — m) arccos (| — l)) 



This completes the proof of the Proposition. □ 
This completes the analysis in the bulk region. We will now move onto the Airy region. 

6.1.3 The Airy region 

The analysis in the Airy region is more difficult compare to the other regions as we will 
need to consider the case where the point t is inside the Airy region. First let us compute 
the asymptotics of the function L n (x) in the Airy region. 

As the asymptotics contain the functions Ai(f n ) and Ai'(f n ), we shall make a change of 
variable and write the asymptotics in terms of the function f n . First, from the expression 
of f n in Definition 16. 1[ we see that the map f n is of the following order inside the Airy 
region. 

U = (£) § ix - 4) (l - ^ + 0((x - 4) 2 )) , (6.38) 
and hence x — 4 is of order 

*- 4 =(9 ,/ -( 1 + s / -© ,+0 (S))- (M9) 

Let us introduce the scaled variable T to be 

T = (j) 3 (t - 4) . (6.40) 

If T is close to the real axis, say T = Tr + iTj, Tr, Tj G M and \Tj\ < 5, then we will 
deform the integration contour R in (L n , L m ) 1 inside a small neighborhood of Tr into the 
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upper/lower half plane depending on the sign of 7j, so that the distance between T and the 
integration contour is always greater than some finite 5 > 0. If T £ R, then for definiteness, 
we will deform the contour into the lower half plane and also deform the branch cut of 
(T — f n )2 into the upper half plane so that the integrals are well-defined. As pointed out 
in the introduction, such deformations will not affect our final result. For simplicity of the 
notations, we will denote the integration contour in the Airy region by R in all cases, but 
with the understand that appropriate deformations of the contours are performed when T 
is close to R. 

As T is of finite distance from the integration contour, we have 



ATj(T-u) T-u V N J 4§JVj(T-ii) 




1+ 4i(A^) + 1 Gi M^h + » g,.,, (6 . 41) 



1 



where Gi and G2 are power series expansions in their arguments with coefficients indepen- 
dent on n and N . Note that the series expansion B starts at power 1 while the expansion 
C starts at power 2. 

Let us write f n = u, then by using (16.4ip and (16.391) . we can write L n (x)/u' as a series 
in terms of u. 

X 1 1 (6.42) 

For some constants do, bo, a k and of the form 

a = 1 + O^" 1 ), 6 = 1 + O (n- 1 ) , a k = a kfl + O (rT 1 ) , 

where a fc and 6^,0 are independent on n. 

As the functions L n are close to each other inside the Airy region, we will introduce 
the following function A(U). 

L N (x(U)) 

A{U) = ^RI7)P (6 - 43) 

where x(U) is the inverse function of U = /W(x). The function A(U) is regarded as a 
function in the variable U. We shall express the double integrals in terms of the function 
A. 
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As in (16.391) we can write x(U) as a power series in U (but with n replaced by N). In 
particular, we have 



where = a>k + O (n -1 ) and bk = bk + O (n" 1 ). If we now let x be the inverse function 
of u — f n (x) and use ( I6.4ip and (I6.42p to represent L n (x(u))/u'(x(u)) as a function in u, 
then we see that, as a function of u, we have 

L n {x{u)) ( Af ^ , Ai{u) (N\~* n M'{u) \ 

\ KJ (T-u)H 4 / (T-u^iVf 2M 7 (6.44) 

x (S n (u)+V n (u)) , 

where Ci >n {u) and V n {u) are power series of the form 

00 p, 

^i,n{U) = / y 2(fc-j+3) ~> Z = ^' 



/ \ \ V dk j i ri u 

V n{u) = 



(6.45) 



2(fc+i)+l , „ . , ' 



where dk,j,i,n an d Ci,k,j,n are bounded and k > 2j. In the series Ci tn (u), both indices and 

j start from 0, while in T> n (u), k and I start from 1 and j from 0. The term S n (u) is of the 
form 

W = 1 - ^ ,n + + O (AT" 1 ) , (6.46) 

where t>o is given by 

43(iV-n) , „, 

v ,n = —± 7 6.47 

m{T-u) 

We are now ready to compute the single integrals. First let us show the following 
Lemma 6.5. Let k,j > 0, then as s± — > +00 and |s_|/|s + | = 0(1), we have 



r ^»= (i.-m-j) +0( i) 

Js_ T — u 2 v / 



(6-48) 



(T-u) 
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Let u(u) be a function of u such that u(u) = 0(u) as u — > oo and let v + = 0(s + ), then we 
have 

/ V 2 / ^dt;dti = s_ tl+fe 2 )+0(l), 6.49 

A_ (T-u)* ■/,,(„) (T-u)f V y 

where Av^> = Ai and Ai^ = Ail . 

Proof. First note that, as the Airy function decays exponentially as u — > oo, 
Ai(u) = -^-=u~^e~^ (l + O («~^)) , u^+oo, 



sin ( !«2 +f 



From this we have 



s+ u k Ai(u 



(6.50) 



u) = — 1 — — fl + O (u i ) ) , u — > -oo. 



-du = 0(1), 



'o (T — u)? 

as s + — > +00. Therefore let us consider the integral in the negative real axis. Integrating 
by parts, we obtain 

u k Ai(u) ^ = f u k j u _ oo A l {v)dv \° _ k r° u k ~ l Ai(v)dv ^ 

s„(T-u^ U { (T-u)i ) s L (T-u)i U (g5i) 

j f° u k Ai(v)dv 
' " x du 



(T-u 



Now by f)6.50p , we see that Ai(w)dw = O u * j as m — ^ —00 and is bounded for u£l. 
This gives us the following estimate 

f°^^du = o(\s^-i-i)+o(i) 

Js- {T-u)i V J 

For integrals involving the derivative Ai'(u), we again note that Ai'(u) is decaying expo- 
nentially as u — > +00 and we again have 



(T — u 



For the integral on the negative real axis, we perform integration by parts again and use 
the estimate 

\Ai(s)\ < C(l + \s\)~* , s<0 (6.52) 
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for some constant C. This shows that the integral on the negative real axis is of order 

f°u k Ai'(u), Ifc _i_i\ 

This proves (16.481) . The estimate (I6.49P now follows immediately from (I6.48P and (I6.50P . □ 

Let us now transform the limits in the Airy region into the variable u. As in the previous 
cases, we are interested in the integral 

I f c+ n fn°+ ~ 

(m n _i) 2 / L n (x)w(x)dx = X 3>n = — / L n (x)dx, 

The following is an immediate consequence of (I6.44p and the estimates ( I6.48p . 

Proposition 6.6. Let T be defined as in \6.J^0^ , then the single integral in the Airy region 
is given by 



iKn-i)' I L n (x)w(x)dx= j M n du + o(N-*\ (6.53) 



N ,. a 
n 



This completes the analysis of the single integral in the Airy region. We will now 
analyze the double integrals. 
By (IBTTOjl . we have 

(x)w(x) / L m (y)w(y)dxdy 



ran 



(6.54) 



= ^ 3 = I- Ln ( x ) / L m (y)dxdy, 

We will change the integration variables to u — f n (x) and v = f m (y). We will denote the 
limits of the outer integral by u± and the upper limit of the inner integral by v + . 

Let us now compute the lower limit of the inner integral. As both u = f n {x) and 
v — fm{y) are conformal inside the Airy region, they can be written as a series of each 
other. Then by (16.391) and the analogue for f m , together with the fact that at the lower 
bound, y = —x, we obtain 

\m / my \n J 20 \n J \n z 

v + — ( — I v 2 + O 

m 



20 V m I V m 2 



The following can easily be seen by writing v as a series expansion in u. 
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Lemma 6.6. Let u = f n {x) and let v(u) be the value of v = f m at y = ^x, then as 
n, m — > oo with m — n finite, we have 



( 4 \ 

u(u) — v + (1 + v i) u + / ] viu 1 , v Q — (n — m)[ — ) + O [m~ 

1=2 W v ' ( 6 .55) 

1 to — n ^ , 9N ^ / 2i+i\ 

<^i = — + (m- 2 , u, = m 3 ), /> 2 , 

15 m v ' \ J 

By (I6.49P and (I6.44p . we see that the double integral is given by 

i^K n -iKm~i / L n (x)w(x) j L m (y)w(y)dxdy = J 3 

inn i r + „ . , N r + „ ., \ i _x\ 

N 3 . 



AT2 



/ + S n A(u) [ + £ m A(v)dvdu )+o( 

Ju- Ju(u) J V 



(6.56) 



Let us consider the terms in f™ + £ n A{u) C\ £ m A{v)dvdu in (16.561) . This term can be 
written as 



/ £ n A(u) / ^ m ^4(f)dvd« 

J u_ J v(u) 

2 , 

= — ~fr ((iV - n) J '(iV - m)\4(j, fc) + (N - n) fc (iV - m) j A(k,j)) ( 6 - 57 ) 

+ o(n-^ s 



for some constants h^j that are independent on N, n and m, where .4(j, fc) are defined as 
follows. 

W> k) = [ U+ P j^%dvdu. (6.58) 

First let us consider the leading order term in (I6.57p . which is given by .4.(0, 0). 
Let C{u) be the following function in u. 

rf n da; JV 2f(-n)s(T-M)^ 

M w ) = i 1 r (6-59) 

d« a*(u)*(4 (t-x*(u)^ 

where x^{u) is defined as the inverse function to u — /jv(^at)- 
Note that is of order 

£W -(t)"'( 1+0 U)) <«°> 
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2 

where the error term is uniform in T as T — > iVX Then A(u) can be written as 

. fn cos rj + (x N )Ai(u) (1 + 0(N~ 1 )) - u~* smr}+(x N )Ai'(u) (1 + O (N- 1 )) 

A{u) - W- -t— -r t~W 

V z (— u) i(l — u) 2 

The following can be derived using (I6.22p . 

2 

Lemma 6.7. The derivatives of A{u) behave as follows as u — >■ —CNz~ £ for some C > 0. 



V2 I (-u)*(T-u)a ' ' V(T-u)3 



(6.61) 



y/2 \ (-u)i(T-u)* \(T-tz)3 
where k > and x^{u) is the inverse of u = /at(x). 

Changing the upper limit t> + into u + , we can write the term .4(0, 0) as 

^ ( °' o) = \ (r + ^ )dM ) 2 ■ § r + ^Tiy^ (z/ - u)j+idM + ° ( e " §uf ) (6 - 62) 

where v is treated as a function of u by using Lemma 16.61 
Let us estimate the order of these terms. 

2 

Lemma 6.8. Let u^/Ns = O (N~ e ). If j is odd, then we have 

J + A(u)A {j \u)(u - uy +1 du = | O (m- 2 +ir) (6.63) 
If j is even, then we have 



4 </«_ (—U)2{T — U) 



+ 



o(iv-f), j = 0. 



(6.64) 



Proof. The first equation can be proven using integration by parts. By (16.551) . (16.611) and 
(16.601) . we have 



A(u)A {j \u)(u - u y +1 du = -A(u_)A [j - l \u_)(u(u_) - u_) j+1 
A'iu)^-^^)^ - uy +1 du + O 



»+ 















(I 




| m~ 2 J 
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From f)6.6ip . (I6.55P and (I6.60p . the first term is of order O [ («_/m3 ] 2 m 2 ) . Repeating 



the integration by parts, we obtain (I6.63p . 

To prove (I6.64p . we have, by (I6.6ip . the following 

A{u)A^\u) - - - — > £ 2 {u) = OU*m >J (6.65) 

as u — y U-. After multiplying by (v(u) — uY +1 and integrate, the error term O (u 1 ^ rn~^ 

gives a contribution of order O (^N~^ for j > and O ^iV -2 ^ for j = 0. 
The term with the cos 2F N factor can be integrated by parts using 

F J , = ^ + + ^ + - i = -(-«). + O^J- i . 
dF w i / 1 

-(-U)a +0 T~T7T 



dw \u2N3 

Then by splitting the interval [w_,0] into [w_,«o] and [it , 0] for some wo between u_ and 
0, and integrating the integral on [it_,iz Q ] by parts, we have 

f \A{u)A«\u) - -C\u) ) (v(u) - uy +1 du 

Ju- \ 4 \ —u) 2 [T — u) / 



max 



(O.j-4) 



O \ j^^m' 2 ] + O (m -2 ^ 



m 3 

Hence we have 

4(-«)l(r-«) 

0fiV-fV j>0; 

o(iv-l), j = 0. 



(6.66) 



As ^(w)./4v)(«)(i/ — ■u}' +1 is integrable in [0, 00) and have exponential decay at +00, f )6.64p 
now follows from (16.661) . □ 

We would now like to combine the terms in f)6.62p with the end point terms J 7 / in ( I6.36P 
from the double integral in the bulk region. 
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Lemma 6.9. Let T n -m be given by \6. 36]) . then we have 



' ir lulu) — u) 



^ L (-u)HT - U)(2k + 1)\ 



k=0 

1 ! n-m)4> If (u{u) -u) 



4 / TT-Ta \7T~ K dx + ~* — C (u)du + O [N~ 



Proof. From (16.361) . we can write T n - m in the following form. (After renaming the inte- 
gration variable from x to xn) 



r4 sin((n-m)0(xjy)) 
Fn-m = 4 / — — tgIxat 



n c _ Nx N (4 — x N )(t — x N ) 



00 H ( -\\kff W W'ik+l ( 6 - 68 ) 



Y> I (-l) fc ((n-m)0(xjy)) 
^ /« Nx N (4 - x N )(t - xat)(2A; + 1)! N 



k=0 



where 4>{xn) = arccos V^- — l). Let u = /jv(^jv)- Then from the proof of Lemma [6.4[ we 
can express 4>(xn) as follows when u — > u_. 

(n - m)4>(x N ) = F„(xat) - F m {^N^j + O (^~f) 
Now F n (x7v) is given by 

7T Il2. .3 7T 

^(zjv) = imp + (x N ) + rj+(x N ) - - = ^gl^) 2 +^7+(^Jv) - ^. 
Similarly, we have F m (xAr) = ~(-/ m )2 + i] + (x N ) - \ and f m {x N n/m) = u(f n ). Then by 

2 2 

Lemma E6] and / n (^jv) = (77) 3 /at(^jv) = (77) 3 we obtain 

_, /n \ 2 /. . ..3 n — N , . ..1 

( n \ 7t ( 1 / u 
+ V+ L Xn )-4 +0 [n U§ 

From 06.81). we obtain 



F n (x N ) - F m (^*n) = (-«) 2 - M ) + f ^ (^j) 

as u — v U-. From f !6.59j> . we see that C 2 {u) is given by 

8(-u)3(T-u) /dx7v N 2 



X^(4 - Xat)2 (t - Xat) 
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du 



By using u = /jv(a;jv), flt>.3j) and (16. ip . we obtain d u x^ = ^-^/xn/xn — 4. Therefore 
C 2 {u) is can be written as 

2 _ 16tt(T - u) dx N 

W ~ Nx N (4 - x N ) (t - x N ) du ' 

This gives us 



4(_1)* ((„ _ m )0( XAr )) 2fc+1 _ 1 u " ( u ( u ) ~ u + ° ( N lu l )) 



2fc+l 

du 



C 2 (u) 



Nx N (4-x N )(t-x N ) 4 (-u)$(T-u) dx^v 

as it — >■ tt_. As v — u = O J and -u_/iV§ = O (iV _£ ), we see that for > 1, 

4(-l) fc ((n-m)^)) 2fc+1 l «*( y («)-tt) W a . ._d«_ / ^ fc -f \ 

iVx A r(4-x JV )(t-x JV ) + 4 (_„)i(r-ti) W da* ^iV 2 ^;' 1 ' 

as m — > w_. Since 4>{xn) = arccos {^f- — l) it behaves as a/4 — a; at as x — )■ 4. Therefore the 
function on the left hand side of (16.691) is integrable for u G (w_,0). Therefore by using 
e < 1/20, we have 



£ c _ iVx(4-x)(t-x) 4J U _ (_ M )^(T- M ) V 

for k > and of order O ^iV _ § J for k = 0. This, together with (16.681) . implies the 
Lemma. □ 

By (I6.62p . (I6.68P and Lemma [6.8} 16. 9[ we see that .4.(0, 0) can be written as 

-y / Mix \ 2 



A(0, 0) - JT n _ m = / A(u)du - / A\u)(u-u)du 

fi u , 7 A " (6-70) 

-4/ » , dx + 0(N-l) 

Jn c _ NxU- x)(t- x) V / 

By (I6.55p . we see that the terms in the sum involving A are of the form 

£ + A 2 (u)(u - u)du = 4 " (? ^T m) £ + A 2 (u)du + O (AH + §) . (6.71) 
Let us now compute the other terms in f)6.57p . These terms are of the form 
((JV - n) j (N - m) k A(j, k) + (N- n) k (N - m) j A(kJ)) 



hkj 



N 3 



for some for some constants hkj that are independent on N, n and m. 
To compute these terms, let us first prove the following. 
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Lemma 6.10. Let F and G be integrable functions on [tt_,it+] such that F and G are of 

3 

order e~ ku ^ as u — > +00 for some k > 0. Then we have 



F(u) / G(v)dvdu + / G{u) / F{v)dvdu 

J y{u) J u— J v(u) 

pu+ r u -\- r u + f u 

= / F{u)du I G{u)du+ I F{u) j G{v)dvdu (6.72) 

J u- J u- J u- J v(u) 

+ [ + G(u) [ F(v)dvdu + (e- ku + j 

Ju- Ju(u) V / 

Proof. The lemma can be proved by integration by parts. We have 

F{u) / G(v)dvdu= / F(u) / G(v)di;du + / F(u) / G(i;)dt;dM (6.73) 

_ J u(u) Ju- Ju J U— J v(u) 

Integrating the first term by parts, we have 

F{u) I G(v)dvdu= / F{u)du I G{u)du 

_ Ju Ju- Ju— 

G{u) J + F{y)dvdu + O (e~ ku + 

o 

as F is of order e~ ku ^ when u — > +oc. This implies 

F(w) / G(t>)dt;d« = / F(u)du / G(u)du - / G(u) / F(t»)dvd« 



y(u) 



+ I + G{u) j F(v)dvdu + O (e 

J u— J v(u] \ 



3 

kul 



'(«) 

This, together with (16.73 p . implies the lemma. □ 

By taking F = A/ (T — -u)- 7 and G = A/ (T — u) k for *4(j, A;), we have the following. 
Corollary 6.1. The terms in ^6.51^ are given by 

(N - m) k ~ 3 A{], k) + {N- n) k ~ 3 A{k,j) = AH{j, k) ({N - mf' 3 - (N - nf' 3 ) 
+ SH(j, k) {(N - m) k -i + (N- n) k -i) 

+ ~((N- m) k ~ j + (N — n) k - j ) J + A 3 du j + A k du + O (e~ ku +^J , 
where AH(j, k) and SH(j, k) and Aj are given by 

Mu) = j0^j]> AH(j, k) = \ (A(j, k) - A(k,j)) , 
SH (j, k) = - ( / Aj(u) / A k (v)dvdu+ / A k {u) I Aj(v)dvdu) 

2 \Ju- Ju(u) Ju- Jv(u) J 
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We can now consider the terms in (I6.57p . We shall first consider the terms where j 
and k are not both zero. From (I6.49p . we see that if j '• + k > 2, then the term will be of 

order O (^N~^. Therefore we shall only consider the cases where j + k < 2. Let us now 

compute the terms A(j, k). 

Lemma 6.11. LetAo(j,k) be 

A (j,k)= / Aj{u) j A k {v)dvdu, (6.75) 

J U— J u 

Then for j + k > 1, we have 

A(j, k) = A (j, k) + 1 ; 1 / Mu) A k (u)du + O (N- 1 ) (6.76) 

N 3 Ju- 

Proof. This can be proved by using the mean value theorem. By repeat use of mean value 
theorem, we have 

/ A k (v)dv = A k {u)(u - v{u)) - Ia'(Q(u - u(u)) 2 (6.77) 

Ju{u) 1 

for some £ u between u and v{u). Therefore we have 
A(j,k) = Ao{j,k)+ / Aj(u)A k (u)(u - u(u))du - - / Aj(u)A' k (Q(u - u(u)) 2 du 

Ju~ * Ju- 

Note that AjA k is of order N^^u^^^ k as u — > — oo and decays exponentially as u — > +oo. 
Then by Lemma [6. 6[ we see that 

/ Aj(u)A k (u)(u-u(u)) i du = —± — AjA k du + 0[N-A , 

Ju- N3 Ju- V ' 

where i = 1,2. Similarly, since j + k > 1, the function Aj(u)A' k (£ u ) is of order N~^u~ 2 as 
u — > — oo and decays exponentially as u — > oo. Hence the integral involving Aj{u)A' k {^ u ) 
is of order AT -1 . This proves the lemma. □ 

From this, ( I6.77P and ( I6.55p . we can compute the terms AH(j, k) and SH(j, k). 

Corollary 6.2. Let AH (i,j) be 

AH (i,j) = -{Mhj)-M3,i)) 
Then the terms AH(j, k) and SH(j, k) are of order 

AH(j, k) = AH (j, k) + (N- 1 ) , SH(j, k) = 1 , ; / A,A k du + O {N' 1 ) 
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By using this, fl6T49|) and ( Egg]) , together with the formula for A(0, 0) (16T70]) and fIBTTTD . 
we arrive at the following. 



Lemma 6.12. Lei n — m = k. Then we have 



J 3 - il n , 3 X mi3 - j;_ m = feJsi + k{2{N -n) + k)J 32 



+ ij n , 3 J A(u)S m du + O (n~ 



(6.78) 



where the terms J 3 j are given by 

4s 4 Z" 4 (/» 4s 
J 31 = r AH (0, 1) r / — — — -dz r / A 2 (u)du, 



4§ / r+ 3 \ 

J32 = 5- / Md« + 7 M 0, 2) 

2iV3 VA_ 4 y 



(6.79) 



As we shall see, the terms J 3 \ and J 3 2 will in fact not enter into the expression of the 
kernel K\. What is important is the structure of equation (I6.78p . 

Remark 6.2. Throughout this section, all the dependence on T are through factors of 
(T — u)~i for j > 0, with the exception of the function C(u). From the definition of C(u) 
in \6. 59\) . we see that it can be written in the form \6. 60\) 

m = l '^\ i (i + o(^ T 



4 J V \Ni 

with an error term that is uniform in T as long as T is not on the integration contour of 

2 

u and v. Therefore all error terms in this section is uniform in T as \T\ — > Ns on the 



contour H + in Theorem 2Jf_ 



6.1.4 The exponential region 

Inside the exponential region [4 + iV _e , oo), we have the following matching formula for the 
polynomials (See Lemma 4.8 of [33J). 

Lemma 6.13. For any d > there exists a constant k > such that, uniformly for 
x G [4 + N~ d , oo), we have 

L n (x) = o(e- fc(:c - 4)n§ ) . 

From this, it is clear that the contribution from the exponential region is of order 
O {e' kN ^ E ^j for some k > 0. 
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6.2 Asymptotics of the skew inner product 

We can now compute the asymptotics of the skew inner products (L n , L m ) 1 . We have 
(L n ,L m ) 1 = - L n (x)w{x)dx L m (y)w(y)dy 

" J — oo J — oo 

POO 

L n (x)w(x) / L m (y)w(y)dydx. 



(6.80) 



By breaking the range of these integrals into different regions, we obtain 

N 2 i^/K n -\K m -i /l \ 1 / 

nm i=i V ^ / ^ v ( 6 .8i) 

for some A; > 0. Then, as X n 2 and X m 2 are of order O ^iV~i j (see (I6.24p ). we see that the 

last term in the above sum is of order O (^N~ t j . By using mean value theorem and the 

asymptotic formulae for A(u), we see that £~ A{u)E m du = O From these, flOTD . 

flCTj) and fl633|) . we obtain 

X n ,3 (X m) 2 + 2m,l) — (X nj l + Zn,2) I m ,3 — X n> 3 / „4.(ll)£ m dll 

J U- 

= Vm^" 1 )^ 3 - (-i) 71 ^) + o (iv-i) . 

From this and (I6.17p . (I6.18p . (I6.37P and Lemma \6.12\ we arrive at the following. 
Proposition 6.7. The product (L n ,L m ) 1 is given by 



nm 



\Lji , Lm) 1 kj 31 + k(2(N - n) + k)J 32 + kJ Br 



Ja + Jat+i (-D-^ + (2(n-W)-*)/ I j, k „M; +Q ^_ S 
l) n + 1 fe/ 1) A; ei>en. 

where are given in \6. 19$ and J B i, J B2 and Ij are given by 



(6.82) 



Io= -ymL A ^ h= -wmhf) L Mu)iu ' 

2tt 2tt 2tt 
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(6.83) 



while J^B3,k is given by 



27r43 T ; et) A 

Jb 3 ,i = 0, J m>2 = — r J-. (6.84) 

We can simplify the expression further by the observation that (L n , L m ) 1 = whenever 
n — m = 1 and n is even. (See Corollary 14. ip Taking k = 1 and n even in ( I6.82p . we have 



- (j 3 i + (2(JV - n) + 1) J32 j + Jbi + Jb 2 + / + (2(n - N) - 1) Ii = O (n~ 

As this holds for any finite N — n as long as n is even, we obtain the following relations. 

J32 = -h + (iV-ft) , -J 3l + J m + I Q = -J B2 + O (iV~i) . (6.85) 

Note that although in the above equation, integration limits u± depend on n, the effect 

of changing these limits will only result in terms of order O (^N~^ J and hence we can 

consider them as fixed under the change of n. 

Let us now compute the factor K n _i/t m _i. We have n n = —-r 2 -. By (14.49 p . we see that 

h n fi/h m ,o — 1 + (N" 1 ). From this, Proposition 16. 71 and (16.851) . we see that the skew inner 
products that we need have the following asymptotics. 

(Lat-i, L N _ 2 ) l = -2/ + 6/1 + O (n~%) , 

(L N , L N - 2 ) 1 = -2I + 2h - J B2 + Jbs.2 + O (n~% 



Remark 6.3. From Ii6.86\) . we see that for large enough N, the products (Ljv-i, £jv-2)i 
and (Lat_3, Ln-a)i will be non-zero if f_° — ^^-du 7^ 0. Since this is an analytic function 

with jump on R that is not identically zero, we can choose the contour E + in Theorem 
such that f 00 — ^-T-dw 7^ onE + . We can therefore assume that f°° — ^-^-du ^ 0. 



Note that, by Remark 16. 2\ these error terms are uniform in T as |T| — > Nz . Let us 
look at the behavior of the skew products when |T| — > and when t remain finite. 

First by (16. 83 p . we see that as T — > 00, the integrals Ij are of the following orders as 
T — > 00. 

/ 2j + l 1 2+j \ 

I- = O (T-—t-zN~—j . (6.87) 

(The statement is clear if \T\ > If \T\ < then by breaking the range of the integral 
into — |T|a) and (— | T\^,u+) and integrate by parts the integral over (w_, — |T|a) using 
(I6.6ip . one can check that (I6.87P is correct.) 
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By using t = 4+ (|) 3 T in (JHSD, we see that Jbz,2 is of the following order. 

J B32 = Q [ N' 1 ( -^-V \ , |T| ATI and t finite. (6.88) 



From these we obtain the behavior of the skew products as T — > oo. 

Lemma 6.14. The skew products (L n ,L m ) 1 in 116.86]) are of the following orders as T — > 



oo. 

2tt 



— ^- (L n , L m ) 1 = O (t-^N-I) , \T\ -»■ Nl and t finite. (6.89) 
hN,o ^ ' 

We can now compute the asymptotic of the kernel. 
6.3 Asymptotics of the kernel Si(x,y) 

First let us used the results in the previous sections to compute the correction term to the 
kernel. Let us define the variables £i and £2 to be 

2 2 

6= (^) 3 (^-4), 6=(^) 3 (y-4). (6.90) 

We will assume that £1 and £2 are bounded from below. First we need to compute the 
integrals e(7r n> iw), which is a linear combinations of the integrals of the L n w. We have 

poo poo 

e{L n w)(y) = - L n (s)w(s)ds + - / L n (s)w(s)ds. 

Jy ^ J — oo 

To compute the first term, let us first assume £2 < N». Then we have, by the asymptotic 
formula of the L n (x) inside the Airy region, and the estimates pQ, 

\Ai(t)\ < Ce-^ t3/ \ \Af(t)\ < C(l +t 1 /4) e -a/3* s/a j t > 0, (6.91) 

the following. 

y/iK n -lN 



n 



poo poo 

I L n (s)w(s)ds= A£rAu + 0[e~ K2 N- 7 A 

Jy Ju(^y) V / 



(6.92) 



for some k > 0, where u = f n {y)- Note that by (I6.9ip . we in fact have an exponential 



decay of order e k ^ in the error term above. However, the decay e k & will be sufficient 
for our purpose. The lower limit u (^y) can be expressed in terms of £2 as follows. 
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where j > 0. By using mean value theorem and f)6.9ip . we can write the integral as 



u { Try) 



2 

1 ( 4 \ 3 



N / 

(N - n) 2 £ 2 (&) + O (e~ k ^/N 7 s 



(6.93) 



2 

for some k > 0. The in (16.931) are given by 

r°° Aif 2 r°° A(v) 

£>=/ Adu+-^A^ 2 ), A = -Afa)- TUTf^du, 

20Ns 2(T-u) 
C - A'(£ ) + M2) + 3 T A{U) du 



(6.94) 



On the other hand, if £2 > N», then from Lemma T6.13I and the asymptotic formula for the 
Airy function, we see that 

poo poo 

/ L n {s)w{s)ds = O (e^ 2 ) , / A£ n du = O (e~ fc?2 ) 



for some constant k > 0. Therefore if we choose the constant k in (I6.93P to be small 
enough, then f!6.93j) remains valid for all £2 bounded below. 
Similarly, the second term in (I6.92p is given by 



L n (s)w(s)ds = ( + ~ j U+ A£ n {u)du + O (N- 

J -00 y2Nt 2 J u _ \ 



Let m_ i jv be /at(c_), then changing the lower limit in the above integral into will only 
result in an error term of order O ^iV~ 1+ ^ = O {^N^i^. Similarly, changing the upper 
limit to +00 will only result in an exponentially small error term. We can therefore change 
the lower limit in the integration to u_,tv an d the upper limit to +00. Let ^ 2 ,n be the value 
of £ 2 at u = w_,jv; then we have 



1 2 



v/m^e (L n w) (y) = * (6) + (jj) % {N - n)*i(&) + \ (J^j " (N - n) 2 * 2 (£ 2 ) ^ ^ 



where ^(£2) is given by 



= - j = 0, 1, 2. 
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Similarly, the orthogonal polynomials L n (x) are given by 

L„(x)w(x) = -A, ( g) * *< 6 ) - g) ' (£±i + („ - *)) 

i \ (6.96) 



where /i and /2 are functions that are independent on N, n and T. Then by using (14.511) . 
we obtain the asymptotics for the correction kernel Ki(x,y). 

y) = (rr — i Qoo + — — — t^Qio 
u -£i) 2 ( T -£i) 2 

jjgOjj^fe) Wfe), A(ft)tt 2 (6.) ^ Agfe)^o(6) n 
+ 1^6)^ Qoi+ (T-fr)i Qn + ^^ Qo2+ (r-fr)i S2 ° (6 ' 9T) 

The coefficients Qjfc are given by 

q 0= (?LY* Q l = ( N \ l N(4h + J B2 - J m , 2 ) | iVT 



4/ 4V2^ ~ V4/ 8V2T/0 \ N J 4V2tF 



.'ivy iV(4/! + Jm-JB3,2) „ fN\~« N 
WiOl — — ~r „ 7= ; , Will 



Q 



02 



4/ 8V27r/ ~' V 4 / 4V2tt' 

/ iV , 3iV(4/i + J B 2 - Jb3,2) 



4 7 V2V2^ 16V2tt/i 



. 4 \ ^ NT fN\~t 3N(4h + J B2 -J B3>2 ) 

^20 — T7 ~ r— + 



(6.98) 



Q 



N{4Ii + J B2 -J B z, 2 ) , /4\3 TAT 



o:-i 



q r - I IV [ _i + T 4 "\ 5 3 ( 4/l + ^ 2 ~ Jm ^ 



N J 8Vt \ 2 \Nj 4/ 
From (16.871) and (I6.88p . we see that as T — > 00, the kernel K\ has the following behavior 

K x {£x, £ 2 ) = O (N%\ , uniformly as |T| JVs and t finite. 
In fact, by using (I6.86p . (I6.89P and (I4.5ip . we obtain the following. 



67 



Lemma 6.15. As T — > oo and t remains bounded, the kernel K 2 and K\ become the Airy 
kernels in the large N limit. 



[jj) a-,(«,,6) = ^ +0 {-T^r 



2 



£2 



lj ifi(6,£ 2 ) = ~^(6) / AH u)<\u + <>(!)< 



(6.99) 



for some k > 0. 

Proof. The statement for the kernel If 2 follows immediately from the representation ( 12. lip 
and the Airy asymptotics of the Laguerre polynomials. To prove the statement for the 
correction term K\, first note that, by ( I6.86p . we see that 

2tt /r _ . 2vr(i-2) 2tt ^ / 25 

(l n , l n - 2 ) 1 = Ar ' -^7 + [N * 



2tT . r , 47T 



47T / 25 \ 

(L N _ 1 ,L N _ 2 ) 1 = — T = = + 0[N-m\ 



uniformly for |T| < ciV3. By dividing the range of integration in f™ + - At ^ 1 du into 
[u-,— \T\%] and [— |T|^,w + ], we see that 

From this, ( I6.100p . (16. 95 p . ( I6.96P and ( 14.5 ip . we see that in this limit, Ki is given by 

2 

3 Kfa, 6) = ^-(6) Q - Ai(«)d« + o (r-f ) + o (at-I)J 

This proves the lemma. □ 



By using the explicit expressions for A and tyj, we obtain the asymptotic formula for 
the kernel when T is finite. 
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Proposition 6.8. Let T be of order Nz c for some positive < c < | ; then for £1 and £2 
bounded from below, the kernel #1(6, £2) is given by 



(i^) 3 Kl ^>& = #1,00(6, £2) + (A^e- fc ^ 

\ «/ — "OO t/ — OO 

Hi / 6 # («)du - # (6) / 2 #i(u)du 



(6.101) 



Ho{Si) 1^ # 2 (u)d W + £ 2 # (6; 



for some k > 0, where Hj, B\ and B 2 are given by 

fc? /J! poo r> poo 1 

B 2 = - -± - - J H du + H,du + - / H 2 d„. 



(6.102) 



As pointed out in [3T], [32] and [33], the eigenvalue statistics will not be affected by the 
rescaling 

#^hr # hr (6-103) 



, 4 / V 4 , 

of the matrix kernel. By rescaling the kernel in this way, all the entries will have the same 
order in the large N limit. From now on, we shall use this rescaled kernel and denote it 
also by K. 

From fl6.10ip . we obtain the following estimate for the rescaled matrix kernel K. 
Proposition 6.9. Let be the 2x2 matrix whose entries are given by 

#00,11(61,6) = #00,22(6,6) = #1,00(6,6) + #2,00(6,6), 

u ft n dK ^ v ft t\ f 6 ur ( t\A (6 ' 104) 

#00,12(6,6) = ^7 , #00,21(6,6) = -/ #11,00 (M,6)dw, 

•96 4 

where #2,00(6,6) is given by 

#2,00(6,6)" 



T-67 6-6 

Suppose T is of order N^~ c for some positive < c < | and that £1 and £2 are bounded 
from below. Let #(6,6) ^ e ^ e rescaled matrix kernel in A6.103\) . then there exists k > 
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such that 



4\i fo(N-Me- k ^) O (N-^e- k ^ 1+ & 



AT / "VSDS^y "OO V"S1 J 'SiJ I / 1 \ 111* 

NJ \ ON' 21) OlN-^e- k t 



Proof. The statement for the 11 th and 22 th entries follows immediately from (I6.10ip . the 
representation (12.111) and the asymptotics of the Laguerre polynomials inside the Airy 
region. The statement for the 12 th entry follows by replacing e(L n w) in (I6.95P by the 
asymptotics of the polynomials in (I6.96P in the derivation of (I6.97p . The computation is 
the same as the derivation of (I6.97P and we shall not carry out the details here. To obtain 
the results for the 21 th entry, we use the fact that e (Si) (x, y) is skew symmetric to obtain 
(See [SI, 122], [23]) 

ry 

e (S x ) (x, y) = e (Si) (x, y) - e (Si) (y, y) = - / S x (t, y)dt. 

J X 

The statement for the 21 th entry then follows from integrating (16.1011) and the asymptotic 
formula for K 2 . □ 

A similar statement can be obtained for the Airy kernels when T — > oo. 

Lemma 6.16. Let K airy (^i,C, 2 ) be the following matrix kernel 

iioin/,11 = K airyM&>£l) = 7 7 ^ o^HSl) / Al(u)du, 

U — ?2 1 J-oo 

d /" 6 

Kairy,l2{£l,&) = —-^rKairy,ll(£,2,£,l), Kairy,21 (£l> £2) = — / Kairy,\\ {U, ^Jdlf. 

Then for T — > oo with t finite and £1 and £2 bounded from below, there exists k > such 
that the rescaled matrix kernel -fr(£i,£ 2 ) in A6.103\) is of the following order as N — )■ 00. 

-J K£ 1 ,&)-K mry fo,S 2 )+^ o(1) o(1)e -fe6 J 

In order to show that the convergence of the Fredholm determinant, we need the fol- 
lowing bounds on the derivatives of the kernel K 2 . 

Lemma 6.17. For T of order N^~ c for some < c < \, we have 



3 

d l d j ( ( 4 \ s 



/or some k > and /, j = 0, 1. 



AT 2 (fi,6) -^2,00(6,6) 
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Proof. The lemma is an immediate consequence of the following results in [32] ([32], (3.8)) 
d l d j ( ( 4 \ s 



2 \ 

(^) 3 K lag (Ci,^) - K 2tairy (Cx,^) J = O (iV-fe-*« 1+ « a )) (6.105) 



where K 2 ^ airy is the Airy kernel 

A 2 ,airy(Cl,6) = ? ? ■ (6.106) 

Si — ?2 

and Ki a g is the Christoffel Darboux kernel of the Laguerre polynomials 

Ki ag (x,y) = w^(x)wl(y)K A N ^ l 



2 L N (x)L N _ 1 (y) - L N (y)L N _ 1 (x) 



x-y 

and x = 4 + £1 (4/N)%, y = 4 + £> (4/iV)^. As K 2 is the conjugate to if Jafl and K 2>00 is the 
conjugate to K 2jairy , 



y{t-y)Y „ „ (T-t 



the lemma follows from (16.105p . □ 

The corresponding statement when T — > oo follows from the same argument but with 
K 2>OQ replaced by K 2Mry . 

With the estimates in Proposition 16.91 and Lemma I6.17[ we can obtain the following 
asymptotic result for the determinant det 2 (I — X^x)- 

Proposition 6.10. Let ( = (z - 4) {N/4}*, & = (a; - 4) (JV/4)s, £ 2 = (j/ - 4) (N/4)*, 



9(0 = a/1 + ^ 2 and G = diag(g,g x ) ; then as N — )■ 00 and T 0/ order up to o ^iV3 



we 



have 



^det 2 (I-XzG(^)KG-^ 2 ) Xz ) = \/det 2 - XcC(6)^ooG- 1 (6)Xc) + 
and for T — > 00 wdw/e t remain finite, we have 



^det 2 (I ~ XzG(^)KG-^ 2 ) Xz ) = v /det 2 (/ - XcG(^)K mry G^(^ 2 )x C ) + o(l), 

where K is the rescaled kernel in \6.10§) and det 2 zs tne regularized 2- determinant det 2 (/ + 
A) = det((J + A)e~ A )e tr( - All+A22 ' ) for the 2 x 2 matrix kernel A with entries Ay. The 
characteristic functions Xz and xc are Xz = X[z,oo) and xc = X[C,°°) respectively. 

The proof of this proposition is exactly the same as the proof of Corollary 1.4 in |32j. 
We shall not repeat the details of the proof here. As in [32J, the function g(£) = a/1 + £ 2 
is to ensure that the 2-determinant exists and there is a great freedom in the choice of the 

<?(£)■ 

We will now analyze the asymptotics of the derivative - lo s ^_ et p in (15. ip . 
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6.4 Asymptotics of the derivative 



d log det 1 
dt 



We will now compute the asymptotics of the derivative dlo ^ et ° i n (j5.ip . Let us write the 
integral in (15. ip in the following form 

9 , , ^ f K 2 (x,x) , f Ki(x,x) , 
— logdetD = -/ A - J dx- K - J dx, (6.107) 



where K2(x,y) is the kernel given by the Laguerre polynomials ( 12. lip and K\(x,y) is the 
correction term on the right hand side of (12.101) . 

To compute the contribution from the kernel K 2 , we will use the following differential 
identity [TB]. (See also [51] Lemma 2.1) Let Z be the matrix related to the matrix Z in 

(14.391) by Z = Zw 2 . Then we have 



Lemma 6.18. Let Z = Zw Q 2 , where Z is the matrix in ^j. 39 ). Then we have 



r K^x)_ dx = i ^-i {t) % (t) \ (6 108) 

J m t — X l V / 



The lemma can be proven in exactly the same way as Lemma 2.1 of The key is to 
use (14.381) . the jump condition of Z on R + and L'Hopital Rule to write K 2 (x, x) as 

K 2 {x,x) = 2~ (tr [zz x Z'_o^ - tr (z+ x Z' + <rty , x e R + 

and then deform the contour of integration to obtain (I6.108p . 

The asymptotics of the matrix Z can be found in [70]. For x ^ [0,4], the asymptotics 
of Z is given by 

Z(x) = (tK N _ l )-^R(x)P OQ (x)e N ^- ip{x) > 3 } (6.109) 

where R(x) is of the form I + O {N" 1 ) and Poo{x) is a matrix bounded in x for x ^ [0, 4]. 
Near the point 4, both the matrix and P^ 1 have a forth root singularity. The derivative 
of R(x) is of order O (iV _1 ) and the derivative of Poo(^) remains bounded. From this, we 
have 



tr (Z- 1 (t) Z' (t) a 3 ) = N fl - yj —j +0(1). (6.110) 

As we will see, this is the part that determines where the saddle point is. Once we have done 
the saddle point analysis later on in this section, we will see that at the phase transition, 
the saddle point will be inside the Airy region. We therefore also need the asymptotics of 
the matrix Z(x) inside the Airy region. This again, can be found in |70j . 
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Lemma 6.19. (170], Section 5.3) The asymptotics of the matrix Z(x) inside a small disc 
Us of radius 5 around 4 is given by 



where rj± is given by Ii6.8\) and A(^) is the matrix 

A(n - v^e-fi ( Az(w2 ° ^ e-f- £ e / 
A ™ ~ VZne \Ai'(0 u 2 Ai>(u 2 0) ' * ' 

A ^- V27re \Ai'^) -Ai'(uO J 6 ' ^ eIV - 

where u = e^ 1 and the regions I, II , III and IV are given by 

J = {£| < arg(0 < 27T/3} , II = {£\ 2vr/3 < arg(0 < vr} , III = 77, IV = 7, 

where the overline indicates complex conjugation. The matrix R(x) is again of the form 
I + O (N^ 1 ). Its derivative is of order R'(x) = O (N^ 1 ). 

From the behavior of the functions r]±, we see that the asymptotics of Z(x) inside the 
Airy region is of the form 



Z(x) = {i KN - 1 )-^R{x)P a {x)N^A{f N {x))w 1 , (6.111) 

where both P a (x) and P~ 1 (x) are bounded and analytic inside Us (See [70J). Moreover, 
as ftf(x) — > oo, there is a matching condition between the formula for Z(x) in the Airy 
region and its formula in the outside region (16.1091) . 

P a (x)N^A(f N (x))w-^ = P^x) (i + O (/; f )) eMf-^H (6.112) 
From this and the fact that is of order O ((x — 4)~^ j as x — > 4, we see that (I6.110p 

2 

remains valid if T — (t — 4) [N/ 4) 3 is large, but with some modifications in the error term, 
tr (z- 1 (t) Z' (t) (7 3 ) =nU- + (n*T~^ + O (nT~^ . (6.113) 

73 



We shall divide the range of T into the regimes < |T| < iV 1 / 5 and \T\ > N 1 ^ 5 and use 
(I6.113P to for |T| > N 1 / 5 . Let us now assume |T| < iV 1 / 5 . Then from (16.1111) we obtain 



tr [Z- 1 (t) Z' (t) a 3 j = tr (A' 1 (f N (t)) A' (f N (t)) f' N a 3 ) 

f 1 -3 , «\ M-N f 2\ (6.114) 
+ tr ^P^P^N^AasA^N-^j +M + O \N~* J , 

The error term is uniform in T throughout £7,5 . By using the asymptotic formula for the 
Airy functions, we see that the matrix A(£) has the following behavior as £ — > oo. 

A(0 = -^r^(_ 1 1 jV^3(7 + 0(r^))e-§«H £^oo (6.115) 
Therefore uniformly in Us, the second term in (16.1141) is of order 

tr [p- l P' a N^Aa 3 A- l N~^ = O ((N/T)*} 
By (16.1151) . the first term in (16.1 14ft has the following behavior as J'n — > oo. 

tr (A- 1 (f N (t)) A' (f N (t)) as) = -2/* + O (/^) , oo. 
Since f N (t) = T (l + O (t/JV* ) Y we have, by mean value theorem, 

tr (A- 1 (f N (t)) A' (f N (f)) /> 3 ) = tr (A -1 (T) A' (T) /> 3 ) + O (t*/N? 

tr (A -1 (T) A' (T)a 3 )+o(rf) 



T 

Therefore after integration, we have 

-t TS ( \ 1 i-T 



[' !^ dx = i /' tr (A -> (T) A' (T) a,) AT + + (N-i 

where |T| < and |T | = For |T| > and |T | = iVs, we have, by (I6.113p . the 
following 



f tr (^(Oi'^aa 



t-4 



7 



= iV / ( I - \/^— I (1/ + O (A'" ) . 

'to ' ' -J to 



Since the first term is of order O (t^j + O ^NzTj, we see that for |T| > the first 
term will dominate over the error term. Summarizing, we have the following. 

2 

Proposition 6.11. Let T = [t — 4) (iV/4) 3 . T/jen t/ie contribution of the kernel K2 to the 
logarithmic derivative of det D is given by 
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1. Uniformly for \T\ < iVs, we have 

f ^^dx = ~ f tr (A- 1 (T) A' (T) a 3 ) dT 

1 x 1 Jto (6.116) 

M / 2 \ 

+ y(f-to) + o(j\rA), 

/or any |To| = , where the integration contour does not cross the jump contours 
of the matrix A. 

2. Uniformly for \T\ > , we have 



7 



for any \T \ = . 

where the integration contour does not cross (— oo,4]. 

Let us now compute the contribution from the correction term Ki(x,y). Before we 
compute the asymptotics, let us make a further simplification using (I4.5ip . First by using 
the recurrence relation of the Laguerre polynomials (14.41 j) . we see that 

Ln = ^IL_^IL = Ln (M _ (tf + jy + i )) _%_^, 
t-x t-x M{t-x) y v " t-x Mt-x 

Substituting this back into (14.511) . we obtain 

Ki(x,y) M ( I Ltv-2 t \ _ (L N ,L N _ 2 ) 1 d 

— dx = -— <- ,L N ) - 2— —— {L N _ 1 ,L N - 2 ) 1 

t-X 2h 0>N -2 \\t-X / 1 (Z/jv-1, i/JV— 2/1 ot 

M d 1 I L N \ N /Ls-i 



h , N _idt 2ho )N -i\t-x / 1 2h 0;N -i\t-x 



(6.118) 



where we have used the fact that (L N , L N ^ 1 ) 1 = 0. 

The main task is to compute the products (L n /(t — x), L m ) 1 . The analysis is very 
similar to those in Section [6j First note that, outside of the Airy region, the analysis in 
Section [6] remains the same and we have 

Lemma 6.20. The single and double integrals in (L n /(t — x), L m ) 1 have the following 
contributions in the Bessel region. 



\lK n —i 



I"' 1 L " (xMx) dx =X ln = ( ~ 1)W f" + O (n-l) . (6.119) 
Jo t-x Jftfi \ J 



^Jnt : i 
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The double integral in the Bessel region is given by 



. , f N 7 LJx)w(x) f N \ . . , . , , 

V^-i^-i / — 7 / L m (y)w(y)dxdy 

Jo t-x J x (6.120) 

(_l)™+n 27r , 5 x 

T/ie single integral in the bulk region is given by 

(«„-,)' f . ^^>Md, = I 2 ,„ = O (n-i) (6.121) 

Let n — m = k, then the asymptotics of the double integrals in the bulk region is given by 
the following. 

j 2 = (n-^), k = 0, J 2 = -^ + (n-^, fc= 1,2, (6.122) 

where Ji is given by 

~ . , f c ~ L n (x)w(x) f c ~ _ . . ... _ 

J2 = V^-i^m-i / 1 — 7 / L m (y)w(y)dxdy. 

Jn~2 t — x J x 



and J 2 is given in {6.31^ 



Note that the choice of e < 1/20 ensures that the error terms are of the same order as 
before. 

In the Airy region, the analysis is more different. Let us now compute these integrals 
inside the Airy region. As in Section 16.1.31 we have 



n 1 r c+ - r v+ / _7\ 

— (iK m -i) 2 / L m {y)dy = / A£ m dv + 0[N b , 

™ J2L X J u{u) V / 

2 



^ ^ r ^ - (t) 1 r ^ + ° 



(6.123) 



where = A/(T — u) and £ n is of the form 

4 = 1- ^o,n + y <„ + O (iV" 1 ) , (6.124) 

where t>o, n is given in (16.47p . By using mean value theorem, (I6.55p . we obtain the following 
estimate. 

v + ( f 4 \ I 



2 



A£ m dv = I A£ m dv + | ( — ) (to - n)A{u)£ m (u) 



(6-125) 

(|) 3 (1 + (I)) + (« - ^) 3 £ m (U 
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for some £ u between u and u(u). As 



A(«)/(y(.-«)t(yd« 



< 



-|T|a 



c 



1(2 (|T| — w) 



-dw 



(1 + |T| 2 )A 



4 | J 

3 



we see that this term is of order O (^N 2 j as T — > 00 and |T| < Applying similar 
argument to the other error terms in (16.1251) . we see that 



A\£ n I A£ m dvdu 

J v(u) 

u + _ r v + ^ / 4 \ I 

A\E n I A£ m dvdu + ( — ) (m — n) J A\AE n E m du 



N 



(6.126) 



4 V (m - n) 2 P+ 



AT 



4 „ 1 
■ 2 



2 

Let us now determine the behavior of these terms as \T\ — > Ns. We have 

Lemma 6.21. As \T\ — > \u-\, the terms Aa(i,j) is of order O ^A^T -1- ^ -1 ^ in T. 

Proof. Let us write the integral in the following form 



A {i,j) 



-It^I /•- |t^| 1 — \t%\ r v+ pu+ pv+ 



Ai 



Ajdvdu 



Ajdu 



Ajdv 

111 j 

T2" 



1 , 
T2" 



Ai 



Ajdvdu 



Then it is clear that the third term is of order O ^N'^T^ 1 ^^ 1 ^ in A and T as \T\ — > \u-\. 
By using the asymptotic formula (16.611) to integrate the second term by parts, we see that 
this term is also of order O ^A^T - * - -^ 1 j . For the first term, we have, from (I6.6ip . the 
following 

i„i 1 i m 1 1 i m 1 

C 



- T J 



Ai 



- IT? 



Ajdvdu 



< 



— \T? 



-du 



N*u(\T\ ~u) i +j+ 1 

for some constant C > independent on T and A. Integrating this gives us the result. □ 
From Lemma 16.211 and the asymptotic behavior of A, we see that the various terms in 

2 

the above expansion are of the following behavior as \T\ — > AX 

f + A^du = O (n-*T-*) , f + A^'du = O (n^T- 2 ) , 

J A\du = O \N-*T~* J , A (i,j) = O\N-3T- i - j - i y 
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From (16.126]) . we obtain the double integral inside of the Airy region. 

Let us now compute the order of the following term in (16.1181) inside the Airy region. 

A= + £w A . (6, 28 ) 

2/i ,jV-2 \ t- X I x 2ft ,JV-l \t-X I x 

First by differentiating the identity (Lat,Lat_i) 1 = 0, we can establish the following. 
Lemma 6.22. Let PH(i,j) be 

i /■«+ /■«+ 

PH{i,j) = - A4u Ajdu-Mhj), (6-129) 
and let J^k be the fallowings. 

J33 = \PH{1, 2) - i (£ + *4?d« - ^ + ^.A'du) , (6.130) 
J M = - (5PH(3, 0) - PH(1, 2)) , J35 = -2 y ^?d«, 

and /e£ Ij be 

71 = ( t J 1 Adu ' /2 = U) vr 1 (6 ' 131) 

Lei k = n — m, then L m } is given by 

( L m \ = -Pff (1, 0) + fcjgi + (iV - n)J 32 + A; 2 jf 33 + (N - nfj M 
h N) o \t-x / 1 

+ k(N - n)J 35 + (-l) n ((-l) fe Ii + (-l) fc (iV - n)h) + J" 2 ( 6 ' 132 ) 



+ O (n~&/ (l + |T|§)) + O (n~*/ (l + Iv^l)) 



The lemma follows from a straightforward calculation using ( I6.126P and Lemma 16.201 
As in the computation of the skew products {L 2 k, L 2k -i) 1 = can be used to simplify 

the expressions (^,L m ) . In this case, we differentiate the identity (L 2 k, L 2 k-i) l = 

with respect to t. Then we obtain 
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By using Lemma I6.22[ we obtain the following identities. 

2.731 - Jz2 - 2I X + 2^ = E u - j u + J 35 -I 2 = E 2 , (6.133) 

at 

where E\ and E 2 are error terms that behave as the error terms in (16.1321) and J 2 in the 
above is taken to be the expression in (I6.37P with n — m = 1. 

We can now compute the term D\ in f)6.128p . By using Lemma 16.221 we obtain the 
following 



^ 2nd t-A d /T . „ fa iia\ 

— = J» - 4J 33 + ?S ^-- S W- + (6.134) 

where the error term E has the same order as the ones in (16.1321) . 
Let us now show that T 2 — 2T\ is of order O (A r_1 ~a). 
From (16.361) , we see that T 2 — 2T\ is given by 

_ „ = 4( 3 m(2^) : 2s,n^) dx = (' 8( S in^(W-l)K 



iv c _ iVx(4 - x) (t - x) ijv c _ iVx(4 - x) (t - x) 

n 

40 3 (1 + (0)) 



-dx. 



'iv c _ 7Vx(4 - x) (t - x) 
where again, = arccos (| — 1J. As <f) — V4 — x (1 + O (x — 4)), we have 



Jiv c _ ivx (t — x) 

n 

2 2 

By the change of variable x = 4 + (]|) 3 £ and t = 4 + (-^) 3 T, it is easy to see that the 
above integral is of order O (iV -1- ^. Hence after integration, we obtain 



D 1 dt = -\ / — v 

2V t 47_ 00 (T- M )2 

3 /- / Ai Ai ^ - (6 ' 135) 



(r-u)*'(T-« 



where (f,g) 1 = L Le(x — y)f(x)g(y)dxdy and the error is of the order indicated in 
(I6.132p . In particular, the error term is of order o(l) uniformly in T = o NV§V For 

T = O [N^j > the error term is of order 0(1). 

Let us now consider the other terms in (16.1181) . Let D 2 be 

M {L N ,L N . 2 \ d M d 

U 2 - — -— ^JV-l, ^7V-2)i - T ~7 {^N+l, L N - 1 ) 1 

ft 0,Af-2 {-L/N-l, <Jt ^0,^-1 Ot 
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By using (16. 86p and ( I4.50p . we see that D 2 is given by 
2itD 2 ( Ah + J B2 - Jb3 )2 , „( AT _*_ m i\\ d 



M 
d 



+ -(4/!- jB 3 ,2 + 0(iV-i)). 



Note that the error term O \ N 24 J is uniform in T as |T| — )■ iVs and its derivative is of 
the same order and is also uniform in T. From this, we obtain 

2nD 2 (A h + J B2 - J B3 , 2 + Q ( N ^\ 8 / + Q / , 
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+ - (4/x - ^3,2 + O (n~ 



By using the definitions of Jo, h and Jbi, Jbz.2 in Proposition 16. 7\ we obtain the following 
for D 2 . 



D 2 dt = J Hi&u - y — log ( t - 3 y # du 

r,oo w ^ \ 2 " " X (6.136) 

-2/ v ~~ 7 dT + 0(iV-AlogT)+0(iV-iT 

2 

The orders of T in the error terms indicates their behavior as \T\ — > iVs. Summarizing, 
we obtain the contribution to the derivative from the correction kernel Ki(x,y). 

Proposition 6.12. Let \T\ < Ns . Then the contribution of Ki(x, y) to the logarithmic 
derivative is given by 



r I. - I .. ..\ I /• n I :2 /• x / /• x 



K\{x x 

_ ' ' dxdt = -7 / ^ l\,^2 dw + / ^ ldM ~ lo S ( / ^odw 





/•°° /4i 2 


4J 


L (T - u) 




77 1 dT- 



oo 



3 /*Y * * ^ ^ r(^^) 2 (6 - 137) 

-- / r, f dT- / J —dT + o(l). 

4 J \{T-u)r (T-u)l) 1 J JZLtfodt* 

where the error term is uniform in T. As T — > oo £mi t remains finite, the function 
P Jr ^t-s^ dxdt u>z// 6e o/ order 0(1) m AT. 

The fact that the integral remains bounded for t finite follows from the estimates (16.127P 
and the expression ( 16.137)) . 
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6.5 Steepest descent analysis 

Before we apply steepest descent analysis to the contour T, let us take a closer look at the 
integration constants in Proposition 16.111 and 16.121 Since the integration constants must 
be the same for all T while the functions involved have jump discontinuities on R, care 
must be taken when choosing these integration constants. 

First note that the expressions in Proposition 16.111 and 16.121 can be simplified further 
through integration by parts. By repeat use of integration by parts, we can write the 
following terms in (16.1351) as 

Ai Ai \ 1 f 00 Ai 2 , 



4 V(T-«)f (T-u)i) 1 2j_ OQ (T-u 



2 J \H X J H 2 dv - ^ '_ - I du + 2 J Hi I H 2 dvdu - \ E x du \ H 2 du 



(Ai'Y-AiAi" (-w)f \ f°° (Ai') 2 - AiAi" 

du — — du. 



T — u 7r(T — u) J J T — u 



i 



The term ~t^~\ is there to ensure the convergence of the respective integrals. From the 
jump discontinuities and the behavior at T — > oo, one can check that 

/ (Ai') 2 - AiAi" f-«)3 \ f°° (Ai'f - AiAi" , 

1 du + / '— du 



T — u tt(T — u) J J T — u 

itr (A- 1 (T) A' (T) as) + 



Therefore we have, for |T| < iVs. 



*-s 7 ./-oo \ J-oo 2n(T-u 

t*T poo pu pT poo poo 

2 I / Hi H 2 dvdudT - / Htdu / H 2 dudT 

Jo J — oo J J — oo J —oo 

/ f'oo jj du \ ^ 

ffidti - log QH H du^j - f I ^ ^ dT + ^ - ?T§ + C + o(l), 

where C is an integration constant that has no jump on R. By Lemma I6.21[ we see that 
the term ( Hidu^j / H du is of order T~f as \T\ > . Hence the we can choose 
the base point of this integral to be +oo. For the other terms, integration by parts shows 
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that 



H 2 dr - , \^ }2 . 1 (1// + 2 / // / tf 2 dnl// 



2vr(T - u) 



'o 




Z^ dM (6.138) 



/o (r-«) 



as T — > oo. Therefore if we let 0± be the origin on the positive and negative sides of 
E, then J °^ <SndT G iR is bounded, where the integration contour goes from + to oo 
along the positive side of R + and then from oo back to CL along the negative side of 
M + . Therefore we can fix the integration constant by defining the integral of <Sn to be 
/ SndT + 1/2 J U+ SndT, which is the same as / S n dT - 1/2 J U+ SndT. We shall use 

these as definitions of the function J T SndT in the upper and lower half planes respectively. 
Summarizing, we obtain 

Proposition 6.13. The logarithm o/detD is given by the following 
1. Uniformly for \T\ < , we have 

SndT - ) 

(6.139) 

x 



det D = C exp ^-S 10 - SndT =f 1/2 jf + SndT - 

/oo 
H du(l + o(l)), ±Im(T)>0 
-oo 

where C$ is an integration constant and Sn is given by A6.138\) . while <S 10 is 

/oo poo I Hidu) o 

H 1 du+ X r ™ UA J dT--T§. 
Jt J^Hodu 3 



2. Uniformly for \T\ > iVs ; detD is given by 




(6.140) 



6.6 Saddle point analysis 



We can now carry out the steepest descent analysis for the contour integral in t. As we 
have seen in the last section, the regularized determinant det2 (/ — xGKG^ 1 ^ remains 
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bounded in N for all values of t. Then from (16.1401) . we see that for \T\ > the saddle 
point is the solution of the following equation. 



2r 



(1 + 

The saddle point is then given by 



0. 



, ._ (1 + r) 

^ saddle 



T 

For r G (—1, 0), we have t sa ddie £ ( _ oo, 0) and for < r < 1/4, we have and t sa ddie G [4, oo). 
Hence the saddle point t sa ddie will not intersect the bulk region (0, 4) for any value of r ^ ±1. 
In this case, we can deform the contour T such that it does not intersect the interval [0, 4] 
and by Proposition 16.10} we see that the kernels Ki and K 2 become the respective Airy 
kernels for all t G T. The largest eigenvalue distribution in this case becomes 

p((A ma:c -4) (AT/4)f <() = C oy J det 2 (/ - XlCoo)GK airy G-^XK,oo)) J (fi(t) + f 2 (t, z)) dt 

for some function /i(t) independent on z and f 2 (t, z) of order o(l). This gives us the Tracy 
Widom distribution for the largest eigenvalue. This is a known result in [39] . 

However, when r = 1, saddle point coincides with the right edge point. In this case, 
the main contribution of the integral in t will come from a neighborhood of t = 4 and the 
kernels will become significantly different from the Airy kernel. This is the case when the 
phase transition happens. Let us now find the steepest descent contour in this case. 

Lemma 6.23. There exists 5 > such that 

<0, |Im(x)|<(5, -<5<Re(x)<4 (6.141) 

Proof. The statement for < Re(x) < 4 is a standard property that follows from the 
Cauchy-Riemann equation. (See, e.g. [T0J ) . Since 





± 



4 - x 



x 



By the Cauchy-Riemann equations, we see that Re ^ (^^7^) ds^ i s decreasing as we 
move away from the real axis. As this real part is zero on [0,4], we obtain the inequality 
(16. 14=1(1 on for < Re(x) < 4. As the function J ^pds behaves as c$ — 2 v / 2(-x)2+0 (x) 
when x — > for some Cq G iR, we see that in a sufficiently small disk D center at x — 0, 
we have 



"X 



Re I / I \/- — - I ds I < °> xeD\ 



1 1 



s 



± 

This completes the proof of the lemma. □ 
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Figure 2: The steepest descent contour T is inside a neighborhood of [0,4] the conditions 
of Lemma 16.231 holds. Inside a neighborhood of 4, we choose the contour to be symmetric 
and denote it by H. 

We can therefore choose our integration contour as in Figure [2] to obtain 

Proposition 6.14. Let w = (f )^ (1 - r) G (-00,00) and let ( = (z - 4) (iV/4)s, then as 
N — > 00, the largest eigenvalue distribution is given by 



for some constant C, where H is a symmetric contour that does not contain any zero of 
J^Hodu and approaches 00 in the sector n/3 < argT < 47r/3, argT 7^ n. It intersects 
the K at a point Tq > £. The function S is given by 



Remark 6.4. In the above formula, the integration should be understood as a sum of 
integrations performed over the contours H + and H_ ; where S± are the intersections o/5 
with the upper/lower half plane. Near the intersection point To, the boundary values of the 
integrand in the upper/lower half planes are to be taken when performing these integrals. 

7 Fredholm determinant 

In this section we will carry out the final part of the analysis and express the determinant 

det2 (/ — X[z,oo)Kx[z,oo)) i n A2.13[) in terms of the Hastings-McLeod solution of the Painleve 

II equation. We will follow the approach in [67]. While the following operations are formal 

and did not take into account the fact that / — GxKooG^x is n °t of trace class, the 

procedure and the resulting formula (17.31) can be justified rigorously as in [BTj . 

Let Si i00 be the operator with the kernel Si i00 = Kx t00 + K 2t oo, D the differential 

operator, e the operator with kernel e(£i — £2) an d X the multiplication by X[c,oo), where 
2 

£ = (z — 4) (A^/4) 3 . Then by Proposition I6.10[ we would like to consider the determinant 




x 



(det 2 (/ - xcG(6)^ooG- 1 (6)xc)) 1 dT, 




5ndT, ±Im(T) > 0. 
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of the following operator. 

j Sl €) xG " (7 - 1) 

where is the transpose of Si )00 and e/ = f(t)dt. Then by (16. 101 [) and (16. 104j) . we 
see that 

d d d d 

- ^-^1,00(6)6) - -^-^1,00(6, £1) = ^--^2,00(^1,6) + ^--^2,00(6,6) 

0^2 C£l C?2 C^l 

= -rff (£i)#o(6) - ^1(6)^1(6) + ^0(6)^2(6) + h 2 (^)h (^). 

This implies Si i00 D = DSf^ and as -De = /, we can follow [67] to write the operator as 

then by using det (J — AB) = det (J — iL4), we obtain 

5 ^-1 1; ^(S 1)00 ) - e X D S* >0O x)[} l) 9 > 

from this, we see that the determinant can be written as the determinant of the scalar 
operator 

det 2 (/ - xGK^G^x) = det 2 (/ - g~ l {e(S 1>00 ) X D + S^x ~ S^xD) 9) ■ 

Let us now show that eSi ;00 = Sf^e. First by writing eSi i00 as 

£51,00 = e (5i )0O ) (Ci, 6) - £ (51,00) (6, 6), 

we see that the first term is equal to Si )00 e because Si jDO D = DSf^, and hence eSi )00 = 
Sf^e. Let us compute the second term. Let / be an L 2 function, then 



/CO POO 
/ 5i )0O (t,6)/(6)d*d6 
-00 J £2 
■t poo POO 

+ ■0 / 5 1)00 (t,e2)/(e 2 )dtde 2 . 

* </-oo ./— 00 

These terms can be computed using integration by parts. The first term becomes 

/oo /»oo /»oo /*^2 

/ s 1 , 00 (jt,&)ffa)dtd&= 5^(6,6)/ /(u)dwifdf 2 
-00 t/ f .1 J— 00 J —00 



(7.2) 



OO /-OO £ 



5i,oo(t,6) / 2 /(u)dudtd£ 2 . 

«/ — 00 
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As - -J-Si,oo(i,6) = §fSi,oo(C2,t), we obtain 



/oo i>co p& roo p& 

/ S lt00 (t,£ 2 ) f(u)dudtd&= 5 1)0O (6,6) / /(«)dudtd£ 2 
-oo »/ £■} J —CO J —CO J — CO 

/co p£ 2 
S , i,oo(6,6) / /(u)dndM6 = 0. 
-00 •/ —00 



-00 »/§2 J —CO J —00 J — 00 

-co f£ 2 



— CO J —CO 



The second term in (17. 2p can be computed similarly and we obtain eS\ :00 = Sf^e. There- 
fore we have 

det 2 (/ - xGK^G-'x) = det 2 (/ - g' 1 {Sfe X D + Six - fixQCD) o) • (7-3) 

As mentioned before, the procedure in obtaining (17. 3p can be rigorously justified as in [67] . 
We shall therefore treat (17.3)) as a rigorous formula and refer the readers to 



7.1 Differential equations for the Fredholm determinant 

Before we compute the asymptotics of this determinant, let us first recall some facts about 
Fredholm determinants and their relations to Painleve equations. 

Let us now recall some basic facts about operators of the form (16. 106[) . The Airy kernel 
(I6.106P is an integrable kernel, that is, it has the form 



K(x,y) = ^ . 7.4 

x-y 

In this case, we have k = 2 and fn = —f 22 = Ai, while f 2 i = /12 = Ai'. Operators of 
these form appear often in random matrix theory and have been studied extensively in the 
literature. These operators were first singled out as a distinguished class in |44J in which 
their properties were also studied. In random matrix theory, they were used to obtain the 
celebrated Tracy- Widom distribution. 

Let us remind ourselves the following well-known facts about integrable operators. (See 
e.g. [30], [44] . |65j ) Let the kernel of an integrable operator on a contour £ be given by 
(17.4)) . Suppose / — K is invertible, then the resolvent R given by 

R = {I-K)- X K= (I-K)- 1 -I, (7.5) 

is also an integrable operator with kernel given by 

„, , SjU Fi j (x)F 2j {y) 

R{x,y) = — , (7.6) 

x-y 

where F^ are given by F^ = (I — K)~ l fij. Let us now specialize to the resolvent of the 
operator with the kernel K 2 ^ ry x acting on R. Then the resolvent of this kernel is given 
by (J — K 2)airy x)~ l = I + RqXi where R has the kernel 

„/,M $0(6)$l(g2)-$0(fr)$l(gl) (7? , 

-ftol4i ? 42j = ~z z • [.'■I) 

?1 — £2 
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where $o and $1 are 

$ = (/ - K^ airyX Y X At, $i = (/- i^^x) -1 Ai'. (7.8) 

For our analysis, we will also need $2 defined by $2 = {I — K2, a iryX)~ 1 Ai". The resolvent 
RqX is closely related to the Hastings-McLeod solution of the Painleve II equation. In fact, 
the determinant det (I + Rqx) is well-known in the random matrix literature and is given 
by the Tracy- Widom distribution for the GUE 



det(I + Rox) = TW 2 ((). (7.9) 

The operator Si )00 has kernel Si i00 = Ki <oc + ip~ 1 K2, a iry4', where ip is the multiplication of 
the function tj) — (T — We will also use the same notation to denote this square root 
function itself. 

By substituting the asymptotic kernels into (17. 3p . we obtain 

det 2 (/ - GxK^xG- 1 ) = det(/ + g- l R T gx) det (/ - g^Kg) , (7.10) 



where K and R are the operators with the following kernels 

K = (I + R T X ) (K^ + K% airy ) (1 - X )e X D + (I + R T x)K1 odX - 



(7.11) 



As the resolvent operator R is conjugate to the resolvent Ro, we have det(/ + g^R gx) — 
det(I + Rox)i which is the Tracy- Widom distribution for the GUE by ( 17. 9p . Let us now 
consider the determinant of I — g~ 1 Kg. We will follow the steps in [H7] to show that the 
operator g~ x Kg is of finite rank, that is, it is of the form 

k 

where / (g) h is the operator with kernel f(x)h(y). As in [ST], the operator (1 — x) e xD is 
given by (See (16) in [BT] ) 

(1 - xhxD = (1 - X) (~e< ® 5 C + eoo <8> , (7.13) 
where e^, €00, 5$ and ^ are given by 

<*(£) = e(£-C), e oo(£) = -~, *c = ^-0, ^oo = ^-oo). 
From these definitions, it is easy to see that 

Cf (l- x ) = Coo (l- x ) = _i(l_ x ). (7.14) 
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Now from (I6.10ip . we see that the kernel K\ )00 is of rank 2. 
#i,oo = H <8> G + H x ® Gi, 
G = (^-I(Ho)-B 1 I(H l )-I(H 2 )+B 2 ), G 1= f^J(#i) + BiJ(ff )^ " 



where J is the integration /(/) = f f(x)dx. By using this, (I7.13P and (I7.14p in (17. lip , 
we see that K can be written as 

K = \{I + R T X ) (ipK^ry^il - X ) + (Ho, (1 - x)) Go 

i 

+ (Hi, (1 - x)) Gx) ® («J C - <W + (/ + # T x) Yl G 3 ® H jX, 



3=0 

where (/, h) = J R f(x)h(x)dx and we have used (a <S> b) (c <S> d) = (b, c)a <S> d to obtain the 
above equation. From this, we see that K is indeed of the form (17.121) with k = 3 and a,, 
/3j given by 

ai = i (/ + R T x) {^K 2Airy r\l ~X) + (Ho, (1 - *)) G + (H u (1 - *)) d) , ^ 
p 1 = 6 ( - 5oo, a,- = (/ + R T x)G j -2, 0j = Hj- 2 X, 3 = 2 > 3 5 

Then the determinant det (i — g^ l Kg s j is given by the determinant of the 3x3 matrix 

with entries 8y — (<Xj, /3j). We will now derive a close system of ODEs in the variable £ for 
these matrix entries. From the expression of Gj in (I7.15p . we see that the matrix entries 
involve the following quantities 

Q_j = ^ (/ + #oX) ^HHj), u-j k = (Q-j, H kX ) , 

V fc = ^(/ + it!oX)#2,a^" 1 (l-x),^x), A; = 0,1, J = 0,1,2. 
We will also introduce some auxiliary functions 

Q+j = (J + RoX) ^I(Hj)i u +>jk = (^ 2 Q+,i, H k x) , 

00 </C 

K K (7.18) 

P+,o = / $ ^d£, P.j = / 

J —DC J —CO 

/CO PCO 
$o^d£, P-j = -J A; = 0,1, j =0,1, 2. 



(7.17) 



which will be more convenient for the purpose of deriving the ODEs. Note that, as in [6 
and [65J, the u±j k can be written as 

u ±jk = ((I + Rox)^ ±l I(H )^H k x) = (^ ±l I(H 3 ),(I + Rox)^H kX ) 

= ^ k ,rl^I(H j ) X ), 9j 

88 



where we have used (17.81) in the above. By using the definition of the resolvent (17.51) to 
write Rqx as R x = K 2 ,ai ry X (I + K^airyXY , we see tnat 



(7.20) 



We have the following relations between these variables. 



Lemma 7.1. The functions Vk and Q+,k in ( |7. i7p can be written as 

V k = V-.k - (H k , (1 - X)) , Q+,j - Q-,j = ~^ u -di - ~^ u -<^ 



R+ ~ R - 



(7.21) 



Proof. By using the definition of the resolvent (I7.5p . we see that (/ + Rox)K2,airy = Ro 
and hence Vk can be written as 

v k = {R^-\i-x)^Hkx) = (4>-\i-x),Ro^H k x) = ((i-x),$k^ 1 )-(H k ,(i-x)), 

where we have used the property of integrable operators (17.81) . The first equation then 
follows immediately from this. 

By using if) = (T — £)a ; one can easily verify the following. 

1 1 $ n $1 $1 $n 

ip tp ip ip 

The equation relating Q±j and 1Z± then follows from this and (I7.19p . □ 

By subtracting 1/2(H , (1 — x)) times the second column, and l/2(ifi, (1 — x)) times 
the third column from the first column of the determinant det (/ — (aj,fa)), we see that 
the determinant is the same as 

A ~lf -("2, A) -fa, fa) \ 
det{I-fa,fa))=det \-\V-, Q I -fa, fa) ~fa, fa) • (7-22) 

\-iv-,i -fa, fa) 1- fa, fa) J 

By using (I7.20p . (I7.15p . (I7.17P and (I7.18p . we see that the other entries fa, fa) can be 
written as 

fa, fa) = \ (l-,o - - Bi ~ ~ (<?-,2 ~ ~ B ^-, 

T 

fa, fa) = — U-fij-2 - BlU-M-l - U_,2j-2 - B2V-J-2, 

f (7.23) 
fa, fa) = ~ (g-,i - + B x ( g _,o - q^f 

fa, fa) = ij-2 + B1U-M-2, j = 2, 3, 
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where q_j are the values of Q-j at ( and q^j are the values of Q-j at oo. 

Let us now derive ODEs for these functions in the variable (. First recall the following 
formula that holds for arbitrary operator K that depends smoothly on a parameter (. 

8 r)K 

-(I - K)- 1 = (/ - K)^ — (I - K) . (7.24) 
Applying this to the operator K 2 ^ airy x, we obtain, as in [65J, the followings. 

0^K 2MryX = -K 2 , mry (^, ()S(C - 6), K 2AiryX y l = -RoiCi, C)P(C, 6) (7.25) 

where p(t;i,t; 2 ) is the kernel of / + R x, that is, p(^i,^) = S(€i - £2) + Ro(£,i, 6)X- Now 
from f)6.106p . we obtain the following. 

This then implies the following for the derivative of p(£,i,!; 2 )- 

Wi + W2) = (/ + ((^ + J2) K ^y^&x) (i + Rox) 



2) 



where we have used the fact that the kernel of the operator (/ <g> g)A is the same as 
the kernel (/ eg) A T g) and that x(I + RoX) is the same as its transpose, together with 

(I + R X)K 2 ,airyX = RqX - 

From this and (17.251) . we obtain the following for p(£,i,£, 2 ) as in [67J. 

r\ r\ r\ \ 

7^7 + 7^ + 7^ p(6> 6) = -*o(£i)*o(&)x, (7.26) 



d( di 2 

From (I7.26P and £2) = <5(£i — £2) + -Ro(£i, 6)X) we obtain the following for £ 2 G (C, °°). 

O O O \ 

+ 77^ + 7^" ^0(6,6) = -$o(£i)*o(6), 6e[C,oo). (7.27) 



d( 0£i «9£ 2 , 

We can now derive a system of ODEs satisfied by the functions in (17.171) and (I7.18p . 

Lemma 7.2. Let q±j and <pj be the values of Q±j and at £ = ( respectively and let the 
function Qj be Qj = q + j — Then functions u±jk, q±,j and Qj satisfy the following 

differential equations. 



Mto*^ * (^V,, ( , 28) 

u-, jk = W(Qj, M- 1 ) /W M- 1 ) , = -ip^q+jfa, 

where W(f,g) is the Wronskian W(f,g) = fg' — gf. 
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Proof. First by using (I7.25P and (I7.19p . one can show as in [65], that, the derivatives of 
the u±jk are given by 



(RfaO^HHJx) <Pk - i^HHj)) (00*, 



From this and (17.211) . we see that 



r> 00 01 

ggj f 00 V f 01 V 



(7.30) 



it. 



jo 



where the prime denotes derivative of (. By eliminating U-ji from the second equation, 
we obtain the differential equation for Qj in (17.281) . From (17.301) . we see that U-jk are 
given by 

u-, jk = W{Q V M' 1 ) /W M' 1 ) ■ 

Now by using the definition of Q+j, we see that 



By applying (I7.26P to this, we obtain 



? +J - + r 1 (C) / ^§^^(6)/(^)(6)d6- 

JR <A 



I M|^ (6m)(6)d6 = _ / ^§^v(6)/(^)(6)d6-0o« + ,io 

= / p(C,6)«7-(^(6i)/(fii)(6))d6-^)«+jo 

JR a t,2 

:Q-,j + 0i ~ °"" ■ jo- 



2^(0 

This, together with (I7.2ip . gives us the differential equation for q + j. □ 

This gives us the first set of differential equations. Let us now derive the second set of 
differential equations for the functions TZ± and Vj. 

Lemma 7.3. Let TZq be 1Z + — 1Z-. Then the functions TZq, 1Z+, V-j and Vqj satisfy the 
following set of differential equations. 

dn o log (0o^- l )^ o _ m w (h$±\ v 



V-j = W (Ko, M' 1 ) /W (0 O ^- 1 , 0i^~ 1 ) 
The functions TZq, 7Z + , V-j and Voj also satisfy the same system of linear ODE. 
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Proof. The proof is similar to ( 17.28]) . First, by (17.25p . we can obtain the derivatives of 
$fc(0 with respect to £. 

^$ fc (O = -i?o(e,O0 fc , 

From this, we obtain the derivatives of P±,o- 



The differential equation for P±,o now follows from the fact that -RoCCijG) is symmetric 
with respect to the interchange of £1 and £2- 

From this and (17.2ip . we again have the following equations for TZq. 



U(o>/ ,x U(0>/ 



This then gives us the Wroskian equations for V-j. Let us now derive the differential equa- 
tion for 1Z + . First by using the fact that Ro is symmetric with respect to the interchange 
of its arguments, we can write R+ as R+ = ip' 1 J Ro(£, C)'0~ 1 (£) c '-£- Then we have 



By using (17.27p . the above becomes 

-K + -ip (0/ V(0 d ^- T77Y^+,o 



ac 2^(0 - - ^ ^ * V(C) 

r 7< -0 T77T' +,0- 



2^(0 u V(C) 

The same argument can be applied to the functions The functions TZq, 1Z + , P-j and Vo,j 
to obtain the same system of linear ODEs. □ 

The ODEs satisfied by the various functions can be further simplified. First, the func- 
tion 0o = {I — K2 ia iryX)Ai(C) is known to be the Hastings-McLeod solution of the Painleve 
II equation (jl.3p (See [66]). 

The Wronskian of 0o/"0 and 4>i/ip can be written as 

W (0 O /V, 0i/V) = (0o0'i - 0i0 o ) = ^- 2 Ro(C, 0- (7-32) 

The function Ro((, C), again is known to be the logarithmic derivative of the Tracy- Widom 
distribution for the GUE [66]. Therefore we have 

/oo 
0o(Od£ = a(C). (7.33) 
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To obtain <pi, let us take the derivative of 4>o and use (17.26p . then we have 



Ai(m = - 



Ai(£)d£-M*o,Aix) 



wi- 



de 

= 0i-0o ($ , Ai X ) ■ 
As in the derivations of (17.29p . we see that 

d d 
—(%,Ai X ) = 4>l —(<!>!, Ai X ) 

Since ($o, Aix) = ($1, A%x) = at ( = oo, we obtain 



01 = 00 + 00 ^ 0o(Od£ = 00 + ^00 

To obtain <p 2 , we use Ai"(£) = £Ai(£) to obtain 

$ 2 (0 = £ (/ + Rox) Ai + $i(0($o, Azx) - $o(0($i, ^x>- 
Therefore we have 

02 



(7.34) 



C + 0o(O0i(Od£) 00-^01- 

Summarizing, we have the following. 
Proposition 7.1. Let {/ be the matrix 



(7.35) 



C7 

















-Wo 



|log(0o/^) 



\-0o 







i 

2i/> 2 



and let hj be the vector hj = ^0, 0, 0, ^-J for j = 0, 1, 2 and hj = for j = 3 and j = 4. 
T/ien i/ie functions in {1.22$ and {1.23$ are given by 



,(oo) 



Hjdu, u. 



q- d ~ qlj = v j4 - v j3 - 

H- = V M - U33 + 1, V-fi = t> 32 , 



,0jr' — U J'2) 

a 

(7 



a 



where Vj is the vector that satisfies the linear system of ODEs 



0,0,0, 



C^+oo, j = 0,1,2, 



dv ■ -> 

-1 = u(C)vj +hj, v r 

v 3 ~ (0,0,-1, Of, C^-oo, tf 4 ~ (0,0,-1, Of, C->+oo 
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These functions can also be characterized using the connection between Fredholm de- 
terminants and Riemann-Hilbert problems. We will outline this connection in Appendix 



Appendix: A proof of the j.p.d.f. formula using Zonal 
polynomials 

We present here a simpler algebraic proof of Theorem 12.11 using Zonal polynomials. Zonal 
polynomials are introduced by James [IB] and Hua [33] independently. They are polyno- 
mials with matrix argument that depend on an index p which is a partition of an integer 
k. The real Zonal polynomials Z P (X) take arguments in symmetric matrices and are ho- 
mogenous polynomials in the eigenvalues of its matrix argument X. We shall not go into 
the details of their definitions, but only state the important properties of these polyno- 
mials that is relevant to our proof. Readers who are interested can refer to the excellent 
references of [55], [50] and [53] . 

Let p be a partition of an integer k and let l(p) be the length of the partition. We 
will use p h k to indicate that p is a partition of k. Let X and Y be N x N symmetric 
matrices and Xi, yi their eigenvalues. Given a partition p = (pi, . . . ,pi( p )) of the integer k, 
we will order the parts p^ such that if % < j, then pi > pj. If we have 2 partitions p and 
p' — (p'ij ■ ■ ■ )Pz(p'))' then we say that p < p' if there exists an index j such that Pi = p\ 
for i < j and pj < p'j. Let the monomial x p be x\ x . . . SEpjS, we say that x p ' is of a higher 
weight than x p Hp' > p. Then the Zonal polynomial Z P (X) is a homogenous polynomial of 
degree k in the eigenvalues Xj with the highest weight term being x p . It has the following 
properties. 



These properties can be found in the references [55], [50] and [M]- Another important 
property is the following generating function formula for the Zonal polynomials, which can 
be found in [50] and pyf] . 



where d p are constants. In particular, if (k) is the partition of k with length 1, that is, 
(k) = (k, 0, . . . , 0), then the constant da.) is given by 



B. 



(tr(X)) fc = 5> p (X), 




(A.l) 




(A.2) 



d(k) 



1 



(2k- 1)!!' 
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For the rank 1 spiked model, let us consider the case where all but one yj is zero and denote 
the non-zero eigenvalue by y. Then from the fact that the highest weight term in Z P (Y) 
is yf 1 . . -yplll], we see that the only non-zero Z P (Y) is Zq^(Y), which by the first equation 
in flA.ip . is simply y k . Therefore the formulae in (IA.1I) and (IA.2I) are greatly simplified in 
this case. 



Jo(N) ik k\ 



x Z (k) (X)y k 

>o(N) ~ k\Z( k )(i N y 

i=l k=0 



(A.3) 



By using the generating function formula, we see that Z^(I N ) is given by 

Z ^ {In) ~ (JV/2- 1)1(2* -1)!!' (A ' 4) 
By taking 9 = ^ in the second equation of flA.3j) . we see that 



(2k-l)\\Z {k) (X)y k 

i=l k=0 



k\ 

fc=U 

We can now compute the integral 



TV 

S(t) 



I e Mt f[(t-x iy )^dt 
Jr i=i 



by taking residue at oo, which is the t 1 coefficient in the following expansion 

i=l k,j=0 J ' 

This coefficient is given by 

N i^Z (k) (X)(2k-l)\\y k M k 



k=0 



Mt -i_ " Z {k) (X)y k M k Mi - / m< XgYgT) 
(I N )k\ (N/2 - 1)\ 



(iV/2-1)!^ Z {k) (I N )k\ (N/2-l)\J 0(N) 

By taking y = 2 (i+ T ) ' ^ ms P roves Theorem 12. II There also exist complex and quarternionic 
Zonal polynomials C P (X) and Q P (X) which satisfy the followings instead. 

11 C P (X)C P (Y) 



f -MReftrfXgYg- 1 )) _ Ml Qp( X )Qp( Y ) 



s P (N) — *! phfe Q p (/jv) 
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Their generating functions are given by 



nc-^r'^|E c ''T' |i ' > 

i,j=l k=0 ' phk v 

n (i - 2^.r 2 =e^e gp(x l Qp(F) 

where C(fc) and g(fc) are 



i,j=l fc=0 phfc ^ p 



2**;!' ^ w (A: + l)!2 fc 

Then by following the same argument as in the real case, we can write down the following 
integral formulae for rank one perturbations of the complex and quarternionic cases. 



[ e-^H^'-^g-^g = / e"< f[(t - x.y)^ it. 



(A.5) 



Appendix B: Connection to Riemann-Hilbert problem 

The functions involved in (I7.18P can also be found by solving a Riemann-Hilbert problem 
(See e.g. [23] . [29] . [30] . [34] . [35J). Let K be an integrable operator acting on R of the form 
(I7.4p and be the column vector with entries /jj. Then operator K\ is also integrable. 
Its resolvent can be found by solving the following Riemann-Hilbert problem. 

1. Y(£) is analytic in C \ [(, oo), 

2. F + (0 =Y.(t)(l-2mf 1 (0/ 2 T (0), ee(C,oo), 

3 . y(0 = / + o(r 1 ), e^oo, 

4. F(0=O(log|e-CI), £^C- 

Then the vectors Fj with entries in (17. 6p are given by i+ (£)/$_,-. 

This Riemann-Hilbert problem can be connected to the Riemann-Hilbert problem of 
the Painleve II equation. We will now outline this connection. For more details, please 
see [29]. Let A(£) be the matrix in Lemma [6.191 Multiplying the solution Y of (1A.1J by 
A on the right and then deform the regions suitably will transform the Riemann-Hilbert 
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problem (lA.lj) into the following Riemann-Hilbert problem. 

1. X(£) is analytic in C \ E, 

2. X + (0 = X_(£)J x , ^S, 



x(0 = 71 r 



(_\ J) e -^(/+o(r 



£ ^ 00, 



(A.2) 



4. X(0=O(log|e-CI), £^C- 
where the contour E is the union of 

S 1 = (-oo,C], S 2 = C + e^M + , S 3 = C + e" 
The contours in E are all pointing towards £. The jump matrix Jx is given by 



J 



x 



1 
1 1 



£ e E 2 u s 3 , J 



o 1 
-1 



Since X = FA, we have the following. 

X n = V^Fe-^(7 - ^.airyX)" 1 ^*, *2i = V^e"^(J - AT^x)" 1 ^'. (A.3) 

Note the difference between our Riemann-Hilbert problem (IA.2[) and the one given in (2.11- 
15) of [29] . In (IA.2j) . as £ — > oo, the next to the leading order term is of order smaller 
than the leading order term. This is due to the asymptotic behavior of the Airy functions in 
the matrix A. As a result, our Riemann-Hilbert problem flA.2[) will be uniquely solvable. 
This is important for us to keep track of the functions $o an d 3>i that appears in the 
resolvent. In [29j, it was shown that the Riemann-Hilbert problem ( 1A.2|) can be solved by 
using the monodromy problem associating with the Painleve II equation. Let E* be the 
union of the contours 



E$,i = e e R + , E^, j2 = e 5 e 1 
and let S±, Si and S2 be the regions 



7T 7T 

S ~ c < ar g£ < 7 
O 



5_ 



77TI 

., Eve, 3 = e e ] 



SvJ/ 



7T 57T 

— < argf < — 
6 &S 6 



5*2 — — Si. 



Then u) is the solution to the following monodromy problem of the Painleve II equa- 
tion. 

1. v) is analytic in C \ E#, 

2. = £eE*, 



4. *(£,v) = { 




1 
±z 1 



£ — )• in S±, ; 

^ in 5i; 
^ in S 2 . 



(A.4) 
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where is given by 

where the contours are all pointing towards infinity. The matrix \I/ also satisfies the Lax 
equation for the Painleve II equation. 

? <9£ V - 2 « - I 4 ^ 3 + (2*9 + W ' 
where q(v) is related to the Hastings-McLeod solution of Painleve II by (see (1.47) of [29] ) 

2-1 (v + 2q 2 {v) + 2q v (v)) = ^u) . (A.6) 

By using the Lax equations (1A.5j) and the behavior of \I/ at £ — )■ oo, one can check that the 
next to the leading term in the expansion in (1A.4j) is given by 

*°° = _1 S (A - 7) 

V 2i 2/ 

The authors in [2S] then showed that the solution to the Riemann-Hilbert problem ( ]A.2j) 
is given by 



(73 



for some h independent on £. The function is determined by making sure that the 
next to the leading order term in X is of order indicated by flA.2[) . By using (|A.7[) , one can 
check that the appropriate choice of h in this case is given by h = —q. The properties of the 
Riemann-Hilbert problem in (1A.4|) is studied very thoroughly in [29], and these properties 
will hopefully be useful in the characterization of the functions in (I7.22p and (I7.23p . 
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